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SECTION  I 
INTRODUCTION 

The  basic  interior  ballistic  problem  is  to  determine  the  energy 
release  and  corresponding  pressure  generated  by  the  burning  propellant 
in  a  variable  volume,  ultimately  to  establish  the  muzzle  velocity  of 
the  projectile.  The  dynamics  must  account  for  certain  losses  which 
include  rotating  band  frictional  effects  and  heat  transfer  from  the  hot 
gaser  to  the  gun.  Krier  and  Adams  [1]  report  that  frictional  losses 
account  for  appr oxitaately  two  percent  of  the  energy  released  by  the 
propellants  in  medium  caliber  guns.  Although  direct  friction  losses 
appear  small,  they  are  important  where  an  accurate  prediction  of 
projectile  velocity  is  desired.  Frict '.on  can  indirectly  influence  the 
thermodynamics  and  heat  transfer  processes  which  play  a  much  larger 
role  in  the  projectile  velocity-  Small  changes  in  the  initial 
engraving  forces  can  increase  the  peak  gas  pressures  and  temperatures  by 
twenty  percent  [2]. 

Another  problem  in  internal  ballistics  is  gun  life.  Both  erosion 
and  wear  can  lead  to  barrel  replaceme^  after  a  few  thousand  rounds 
[3]-  In  the  past  rotating  bands  have  been  made  from  bronze  (90  CU,  10 
ZN),  called  gilding  metal;  however,  recent  success  [4]  with  nylo»  bands 
suggest  that  they  will  be  used  extensively  in  the  future.  .  Plastic 


bands  not  only  reduce  sliding 
for  b 


forces  but  Improve  gun  life.  ■  The  need 
er'  velocity  predictivfe  codes  and  for  an  understanding  of  barrel 
wear  has  intensified  Interest  in  rotating  band  contact  loads  and 
sliding  forces. 

Frictional  behavior  in  internal  ballistics  is  exceedingly  complex 
due  to  the  large  loading  forces,  high  sliding  velocities,  and  the 
nature  of  the  dynamically  changing  interface  between  the  projectile  and 
barrel.  Engraving  and  bore  sliding  force  equations  take  essentially 
two  forms.  In  the  closed  form  models  [5]  the  friction  is  assumed  to  be 
proportional  to  the  kinetic  velocity  of  the  projectile.  In  open  form 
models  a  table  of  resistance  force  as  a  function  of  projectile  position 
is  used  or,  more  simply,  a  corstant  engraving  force  followed  by  a 
smaller  constant  bore  force  [(].  Estimates  of  these  two  forces  are 


given  as  10J  and  if  of  the  maximum  gas  pressure.  Recently,  Fisher  and 
Trippe  [7]  have  divided  the  friction  force  into  a  linearly  increasing 
force  during  the  engraving  process  followed  by  a  linearly  decreasing 
force  during  the  bore  sliding  process.  This  model  is  based  upon  data 
obtained  from  extruding  brass  and  aluminum  stock.  In  all  cases 
empirical  sliding  forces  are  based  on  ad  hoc  conditions. 

1.  THE  PROBLEM 

Recently  [83,  the  writer  has  developed  a  theory  fcr  the  normal 
loads  and  sliding  forces  encountered  by  projectiles  when  their 
rotating  bands  (plastic)  are  engraved.  The  theory  illustrated  the  need 
to  consider  geometrical  details  and  served  to  isolate  the  dynamic  flow 
pressure  and  the  coefficient  of  friction  as  the  main  parameters.  The 
theory  also  accounted  for  band  radial  displacements.  These 
displacements  are  caused  by  normal  loads  on  flexible  walled  projectiles 
and  by  propellant  gas  pressures  which  produce  radial  strain  in  the 
projectile  and  barrel.  Finally,  the  theory  was  extended  to  determine 
the  normal  loads  between  the  projectile  and  barrel  for  thr 
pC3t-engraving  region  of  contact.. 

An  application  of  the  theory  requires  a  knowledge  of  sliding 
friction  as  a  function  of  load  and  velocity.  Generally,  this  ,  . 
information  is  available  only  for  selected  materials  and  loads.  The 
effect  of  velocity  on  the  coefficient  of  friction  has  been  one  of  the 
more  neglected  areas  of  research  in  tribology. .  In  most  cases  the 
coefficient  of  friction  decreases  with  velocity.  Above  a  certain  PV 
(pressure  x  velocity)  the  coefficient  of  friction  drops  to  a  value 
substantially  less  than  0.1  and  remains  constant.  The  critical  PV 
value  for  bronze  rotating  bands  is  approximately  4  x  10**  (pslHfps), 
and  Montgomery  [9][10][11]  identifies  this  equilibrium  region  with 
material  melting.  Since  polymers  melt  at  much  lower  tea^eratures  than 
bronze,  it  is  probable  that  most  of  the  band  engraving  takes  place  in 
this  equilibrium  region.  When  Stiffler's  engraving  theory  £8]  was 
applied  to  the  experimental  work  of  Cross,  it  was  found  that  the 
coefficient  of  friction  for  plastic  bands  was  on  the  order  of  0.02. 

It  is  proposed  to  develop  a  theory  of  friction  for  plastic 
rotating  bands  in  interior  ballistics.  Section  II  reviews  the  present 
understanding  of  friction  in  the  scientific  community.  Section  III 


reviews  the  rheology  of'  both  Newtonian  and  non-Newtonian  fluid 
mechanics  since  the  underlying  phenomenon  centers  about  a  melt  layer 
between  surfaces.  Section  IV  develops  a  new  melt  friction  theory  which 
is  applied  to  metal  rotating  bands  (Newtonian  melt).  The  theory  is 
used  to  explain  the  rotating  band  data  of  Montgomery  [9-10].  This 
section  is  a  self-contained  paper  [12]  accepted  for  publication  in  the 
ASME  Journal  of  Tribology.  Section  V  extends  the  theory  to  the 
circular  contact  of  pin-on-disk  devices  which  are  the  main  source  of 
friction  data.  A  comparison  with  the  high  speed  Franklin  Institute 
data  [11]  further  validates  the  theoretical  approach  go  plastic 
rotating  bands.  This  chapter  has  been  submitted  for  review  to  the 
Journal  of  Tribology  [131.  Section  VI  develops  the  theory  for  plastic 
rotating  bands  which  requires  a  major  adaption  for  non-Newtonian 
polymer  melts. 
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SECTION  II 

LITERATURE  REVIEW  OF  FRICTION 

).,  FRICTION  (GENERAL) 

The  basic  mechanism  of  friction  is  non  well  established  [14], 
particularly  for  metals.  It  is  known  that,  when  two  surfaces  are 
placed  together,  intimate  contact  occurs  only  at  the  tips  of  the 
asperities  in  the  surface.  These  asperities  are  then  deformed 
plastically  until  the  real  area  of  contact -A*,  is  sufficient  to  support 
the  load  W.  Over  this  contact  area  there  is  marked  molecular 
interaction  so  that  a  cold  weld  or  junction  is  formed  which  is 
comparable  in  strength  to  the  bulk  material  forming  the  junctions.  In 
sliding,  a  certain  force  F  is  required  to  shear  these  minute  junctions. 
Thus, 

F  -  Ars  (1) 

where  s  is  the  shear  3trength  of  the  Junction.  Since 


where  Pg,  is  the  material  flow  pressure,  the  coefficient  of  friction 


Although  this  model  explains  the  main  characteristics  of  metallic 
friction,  the  value  of  the  friction  coefficient  calculated  from 
material  strength  properties  are  somewhat  lower  (4  ar  0.2)  than  the 
value  found  in  practice  (4  s  1.0).  The  reason  for  the  discrepancy  was 
first  explained  by  McFarlane  and  Tabor  [15].  As  sliding  commences,  at 
an  asperity  Junction,  the  real  area  of  contact  can  increase  several 
fold.  This  is  because  the  condition  for  plastic  yielding  of  a  Junction 
is  determined  by  the  combined  effect  ot  the  normal  stress  p  and  the 
tangential  stress  3,  the  yield  criterion  being 

p2  +  -^2  .  pQ2  (4) 
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where  p0  is  the  initial  static  contact  pressure,  and  Y  is  a  constant 
with  a  value  of  about  ten.  As  soon  as  a  tangential  force  is  applied,  p 
diminishes  with  a  corresponding  increase  in  the  contact  area. 

2.  POLYMER  FRICTION  (LOAD) 

One  of  the  first  studies  of  polymer  friction.  Shooter  and  Tabor 

[16] ,  showed  that  strong  adhesion  occurs  at  the  sliding  surface  and 
that  over  a  restricted  range  of  loads  frictional  behavior  is 
approximately  expressed  in  terms  of  equation  (3).  Howell  and  Mazur 

[17]  pressed  polymer  fibers  against  a  flat  plate  and  reported  that 
friction  does  not  increase  linearly  with  load: 

F  -  Wn  (5) 

where  n  <  1 Thus,  the  coefficient  of  friction  decreased  a3  the  load 
increased.  They  interpreted  this  result  as  meaning  that  the  real  area 
of  contact  followed  a  law  of  the  type:  A  ~  Wn.  Pascoe  and  Tabor  [18] 
investigated  the  friction  of  crossed  polymer  fibers  (diameter  D  mm)  at 
very  light  loads  and  confirmed  that 

u  -  W-m  (6) 

(Figure  1).  In  contrast  to  metals  there  appeared  to  be  little  or  no 
junction  growth  during  sliding.  No  distinction  was  made  between  the 
real  contact  area  and  the  apparent  contact  area.  A  very  thorough  study 
of  nylon  friction  was  carried  out  by  Adams  [19][20].  He  slid  nylon 
6-10  hemispheres  (diameters:  0.24  -  1.16  cm)  on  a  smooth  glass 
surface,  using  speeds  from  10”7  to  3  x  10-2  cra/sec  and  loads  from  0.7 
to  200  g.  His  measurements  showed  that  the  load  dependence  of  both  the 
apparent  contact  area  and  the  friction  force  could  be  represented  by 
the  expressions: 

A  -  0W®  (7) 

and 

F  »  aWn  (8) 

where  8  and  a  are  constants.  The  values  of  m  and  n  were  independent  of 
specimen  radius  and  were  given  by  0.708  and  0.781,  respectively  Adams 
proposed  that  the  increase  of  F/A  with  load  could  be  explained  either 
by  an  Increase  in  the  ratio  of  true  to  apparent  contact  area  with  load 
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Figure  1.  Friction  of  Fibers  as  a  Function  of  Load 


while  the  shear  strength  of  the  true  contact  area  remained  constant,  or 
by  an  increase  of  this  shear  strength  with  pressure.  After  detailed 
considerations,  he  concluded  the  latter  was  more  likely. 

3.  AREA  OF  CONTACT 

Ic  was  initially  accepted  that  all  contacting  asperities  deform 
plastically  since  Amontons'  laws  of  friction  woula  be  violated  if  the 
contacts  were  elastic,  i.e.,  Ar  -  W2/3  from  Hertzian  contact  of  a 
single  asperity.  However,  Lodge  and  Howell  [21]  and  Archard  [22] 
proved  that  multiple,  contacts  with  elastic  asperities  will  produce  a 
real  contact  area  which  is  nearly  proportional  to  the  load.  Lodge  and 
Howell  assumed  a  spherical  surface  of  very  small  hemispheres  of 
constant  height  (Figure  2).  For  elastic  deformation  the  real  area  of 
contact  turned  out  to  be  Ar  -  Archard  allowed  for  asperities  of 

variable  heights  and  found  that  A  -  W*®&.  Furthermore,  Archard  showed 
tnat,  if  the  asperities  are  covered  by  smaller  asperities,  A  -  W2^/27. 
In  general,  both  elastic  and  plastic  deformations  probably  occur. 

Plastic  deformations  predominate  with  rough  surfaces  and  large  loads. 

For  polymers  subjected  to  small  loads  the  apparent  contact  area 
between  a  hard  surface  and  polymer  sphere  is  given  by  A  -  w(2m-m  ). 
Thus,  real  areas  of  contact  for  polymers  are  essentially  proportional 
to  the  load  even  when  viscoelastic  deformation  of  asperities  occur. 

4.  POLYMER  FRICTION  (LOW  SPEED) 

In  the  conventional  presentation  of  the  adhesion  theory  of 
friction,  it  is  assumed  that  the  small  contact  area  increases  with  the 
load  in  a  way  that  the  average  contact  stress  remains  constant  and 
equal  to  the  flow  strength  of  the  softer  material,  then  the  3hear 
strength  of  the  junctions  can  be  treated  as  Invariant.  Thus,  the 
friction  coefficient  is  constant  and  independent  of  load  as  required  by 
Amontons'  law.  However,  polymers  do  not  obey  Amontons'  law.  A  number 
of  different  workers  [24][25][26][27][28] ,  taking  the  lead  of  Adams, 
verified  that  the  strength  properties  of  polymers,  unlike  metal3, 
increased  with  the  hydrostatic  pressure.  Towle  [29][30]  proposed  a 
simple  extension  of  adhesion  theory.  If  the  shear  stress 

s  -  s0  ♦  op  ,  (9) 

it  follows  that 

'  ’  8 


Figure  2.  Model  of  Surface  Roughness 


s  *  sQ  +  ap  . 
it  f ol low3  that 

y  »  —  ♦  a  (10) 

P 

The  equation  predicts  that  the  friction  coefficient  decreases  with 
increasing  load  (pressure),  approaching  an  asymptotic  value  at  high 
loads.  Towle  C30]  states:  "for  many  metals  the  parameter  a  is  very 
small,  typically  7O.OI,  while  for  polymers  it  is  sometimes  an  order  of 
magnitude  larger.”  He  found  that  a  «  0.07  for  the  phosphinate  (Figure 
3).  The  agreement  between  experimental  values  of  friction  coefficient 
and  the  above  model  is  exdellent.  At  low  loads  the  area  contacts  are 
partial,  but  at  high  loads  the  real  area  of  contact  approaches  the 
apparent  area  of  contact.  The  contact  stress  pc  at  which  the 
transition  from  low  load  to  high  load  occurs  can  be  shown  to  be  [30] 


According  to  plasticity  theory  the  bulk  flow  pressure  pm  should  be 
about  five  times  the  bulk  shear  strength  of  the  material.'  Then  the 
transition  stress  i3 


Petersop  and  Ling  [31 ]  suggested  a  relationship  similar  to 
equation  (10)  for  metals.  Generally,  the  coefficient  of  friction  was 
independent  of  load,  obeying  Amontons'  law,  until  loads  were 
sufficiently  high  to  produce  gross  deformations  as  with  metalworking 
processes.  Their  values  of  a  -  0.11  for  aluminum  contradicts  the 
statement  of  Towle  that  a  is  small  for  metals.  Undoubtedly,  the 
transition  stress  for  metals  occurs  at  levels  outside  of  typical 
engineering  applications  and  experimental  friction  studies. 

Briscoe  and  Tabor  [27]  studied  the  effect  of  pressure  on  the  Shear 
properties  of  very  thin  films.  They  found  that  the  shear  strength  of 
thin  films  is  much  less  than  the  shear  strength  of  the  bulk  polymer. 
However,  the  asymptotic  friction  coefficient  a  showed  little 
difference.  It  is  apparent  that  the  constants  s0  and  a  can  vary  with 
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Figure  3.  Shear  Strength  and  Friction  Coefficient 
of  Poly  (Chromium  Phosphinate) 


sliding  velocity  since  contact  temperatures  increase  with  frictional 
energy  dissipation.  Briscoe  and  Tabor  [27]  reported  that  the  shear 
strength  of  polymethylmethacrylate  (PMMA)  abruptly  dropped  when  test 
temperatures  reached  the  glass  transition  temperature  (Figure  4). 

Ludema  and  Tabor  [32]  found  a  pronounced  variation  in  the  sliding 
friction  of  polymers  with  temperature  and  speed  (less  than  1  cm/sec). 
There  is  a  dramatic  decrease  in  junction  shear  strength  s  and  increase 
in  contact  area  A  for  temperatures  up  to  150°C  (Tigure  5).  The  product 
sA ,  which  is  the  friction  force,  leads  to  unpredictable  behavior  in  the 
coefficient  of  friction.  The  effect  of  speed  (10-1*  to  1  cm/sec)  on  the 
coefficient  of  friction  is  minimal  for  nylon  6-6  but  13  substantial  for 
polythene  (Figure  6).  The  reader  may  contrast  the.  above  speed  effect 
on  nylon  6-6  with  the  results  of  Adams  for  nylon  6-10  (Figure  7). 

5.  POLYMER  FRICTION  (MEDIUM  SPEED) 

When  two  surfaces  are  in  relative  motion,  the  frictional  work  will 
be  liberated  as  heat  at  the  surface  asperities.  Thus,  temperature 
effects  should  be  more  pronounced  at  higher  sliding  speeds.  Bowden  and 
Tabor  [14]  obtained  experimental  evidence  for  surface  contact 
temperatures  as  high  as  1000°C  with  constantan  sliding  on  steel.  With 
other  metals  the  peak  temperatures  recorded  corresponded  to  their 
melting  points.  It  was  noted  that  .lubricants,  although  reducing  the 
peak  temperatures,  did  not  prevent  hot  spots  from  occurring.  Also,  the 
sliding  of  poor  conductors  resulted  in  correspondingly  higher  surface 
temperatures. 

From  elementary  considerations  the  rate  of  friction  work  must  be 
equal,  to  the  rate  of  heat  transfer  to  the  material  bulk; 

WUU  -  cA(AT)  (13) 


where  c  is  a  constant  which  is  a  function  of  the  thermal  conductivity 
k.  I'  the  temperature  change  is  limited  to  the  material  melting  point, 
the  above  equation  car  be  rearranged  to  the  form;  pressure  x  velocity 
-  constant,  the  limit  equation  for  bearing  failure. 

Jaeger  [33]  has  given  an  equation  for  calculating  the  surface 
temperature  of  a  single  asperity ,  contact  of  radius' R0; 


AT 
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Shear  Strength  as  a  Function  of  Temperature 
for  PMMA  Films  at  a  Contact  Pressure  .of 
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Coefficient  of  Friction,  Area  of  Contact 
and  Shear  Strength  for  a  Hard  Steel 
Slider  on  (POM)  and  (PP)  as  a  Function 
of  a  Temperature.  Sliding  Speed  0.035  cm/ sec 


Following  Jaeger,  Archard  [3*0  found  that,  for  plastic  deformation  and 
low  speeds  (L  <  . 1 ) , 


and  that,  for  plastic  deformation  and  high  speeds  (L  >  100), 

,1T  .  U’/g  ■  „ 

3.25(pck)1/'2 


where  L  »  R0U/2a;  a  *  thermal  diffusivity;  c  »  specific  heat. 

Vendgradov  et,al.  [35]  made  a  study  of  thermal  effects  by  sliding 
polypropylene  disks  on  metal  rings.  The  conditions  were  described  as 
heavy  duty  although  the  sliding  velocity  was  only  14  cm/se,c.  Figure  8 
shows  the  coefficient  of  friction  for  dry  polypropylene  a3  a  function 
of  load.  The  curves  have  a  maximum  which  was  identified  with  surface 
melting  of  the  polymer.  Increasing  load  at  constant  sliding  velocity 
caused  an  increase  in  friction.  Thi3  was  attributed  to  an  increase  in 
the  contact  area.  At  a  critical  load  the  polymer  begins  to  melt,  and 
the  friction  coefficient  decreased  with  the  increasing  load.  A  thin 
melt  layer  was  foi  id  on  the  surface'  for  the  decending  branch  of  each 
curve.  This  layer  was  detected  by  rapidly  cooling  the  sample  and 
observing  the  newly  formed  hard  amorphous  layer.  Test  results  were 
grouped  into  two  curves  according  to  the  different  thermal 
conductivities  of  copper  and  aluminum  on  the  one  hand  and  of  steel  and 
cast  iron  on  the  other. 

.McLaren  and  Tabor  [36]  presented  one  of  the  few  studies  of  polymer 
friction  as  a  function  of  typical  engineering  speed  (Figure  9).  The 
results  show  the  distinct  maximum  coefficient  of  friction  which  occurs' 
at  a  given  speed.  Although  they  speculate  that  the  temperature  rise 
was  probably  of  the  order  of  tens  of  degrees  C  at  less  than  100  em/sec, 
it  is  certainly  greater  at  higher  speeds. 

The  friction  of  nylon  against  steel  was  investigated  by  Cerico 
[37]  with  an  apparatus  that  used  rotating  disks.  As  the  disks  rotated 
-;t  different  speeds,  a  sliding  speed  of  76  cra/sec  was  produced  in  the 
contact  zone.  The  coefficient  of  friction  required  a,  run-in  time 
before  it  stabilized  (Figure  10).  The  average  roughness  depressed  with 
time  until  it  reached  30f  of  its  initial  value.  As  others  have 
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Figure  10.  Coefficient  of  Friction,  Ovalization, 
Temperature  Rise,  Wear  of  Nylon  AQ  and 
Total  Wear  AQ'  of  Nylon  on  Steel  Pairs 
vs.  Average  Stress  in  Contact  Area 
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observed,  an  increase  in  the  friction  coefficient  occurred  up  to  a 
maximum  beyond  which  a  noticafcle  decrease  took  place.  The  peak 
friction  coefficient  was  linked  to  a  significant  rise  in  measured 
surface  temperatures,  probably  well  below  the  asperity  contact 
temperatures.  The  friction  coefficient  decreased  with  the  formation 
of  a  layer  of  nylon  on  the  steel  surface.  This  layer  melted  for  the 
higher  loads.  Just  recently  Clerlco  [38]  reported  on  friction  and  wear 
studies  of  polyacetals  at  a  sliding  speed  of  0.45  m/sec.  The  friction 
coefficients  again  displayed  a  maximum  at  a  load  of  50  Newtons  and 
approached  an  asymptotic  value  of  0.1  (Figure  11).  Wear  of  the  much 
harder  glass  filled  polyacetal  (GV  1/30)  wa3  at  least  40  times  tne 
unfilled  material,  counter  to  wear  theory.  As  observed  earlier  with 
nylon  composites  and  non-composites  [37],  the  surfaces  appeared  similar  ' 
to  desert  sands  which  confirmed  the  presence  of  polymeric  melt.  The 
melted  layer  solidified  into  an  amorphous  brittle  layer  with  a  softer 
sublayer.  The  thickness  of  the  brittle  layer  depended  on  the  normal 
load  and  sliding  distance.  For  non-composites  the  layer  was  about  3-5 
um  thick  at  W  -  50  Newtons  and  about  15  um  thick  at  W  *  200  Newtons. 
Catastropic  wear  at  high  load3  and  long  sliding  distances  was 
identified'  with  the  disruption  of  this  brittle  solidified  layer. 

Tanaka  [39][40]  investigated  the  friction  and  wear  of  both  glass 
and  carbon  filled  polyacetal  and  teflon  at  sliding  speeds  up  to  2.5 
m/sec.  Polymer  pins  (3  mm  dla.)  were  rubbed  on  both  steel  and  glass 
disk3.  It  was  found  that  melting  of  the  frictional  surface  layer 
occurred  easily  under  ordinary  rubbing  conditions..  The  melting  depth 
of  the  unfilled  polyacetal,  rutbing  against  steel,  was  generally 
several  microns  thick  while,  ,'n  the  case  of  a  glass  disk  at  higher 
speeds,  it  was  about  twenty  microns.  For  filled  polyacetals  the 
melted  layer  was  associated  with  a  very  dense  fiber  rich  surface.  The 
carbon  filler  produced  a  matted  layer  over  50  um  thick  (Figure  12) 
while  the  glass  filler  formed  a  mat  only  near  the  center  of  the.  pin 
specimen  which  was  attributed  to  the  poorer  thermal  conductivity  of 
glass. 

The  coefficient  of  friction  of  the  filled  polymers  against  glass 
and  steel  is  shown  in  Figure  13-  The  friction  coefficient  was 
considerably  higher  than  the  values  reported  by  Clerlco  [38]  for 
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Figure  11.  Wear  Volume  and  Friction  Coefficient  of 
Plastics  vs.  Normal  Load 
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Figure  12.  Scanning  Electron  Micrograph  of  Carbon 
Fiber-filled  Polyacetal  after  Rubbing 
under  a  Load  of  50  N  against  a  Glass 
Disk  at  2.5  m/s 
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Figure  13. 
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identical  loads  and  speeds.  The  difference  apparently  was  due  to  the 
contact  stress  (recall  section  3):  Tanaka  (7  N/mm^) ;  Clerico  (11  to 
65);  unfilled  polyacetal  yield  strength  (73)-  The  fillers  promoted 
lower  friction  because  o^  the  reduced  adhesion  between  fibers  and 
steel.  The  differences  were  less  pronounced  for  the  glass  disk  since 
glass  limited  heat  conduction  to  the  disk  and  caused  melting  at  lower 
speeds.  In  the  non-melt  region  below  0.5  m/sec,  the  unfilled 
polyacetal  had  a  much  lower  wear  rate  than  the  filled  polyacetal.  At 
higher  speeds  the  filled  polyacetals  showed  marked  superiority.  Tanaka 
attributed  this  behavior  to  the  ability  of  the  fibers  to  prevent  an 
outflow  of  molten  materials.  However,  it  is  also  likely  that  any 
solidified  fiber-filled  layer  has  a  much  higher  strength:  than  the 
friable,  pure,  polyacetal  layer. 

6.  POLYMER  FRICTION  (HIGH  SPEED) 

Friction  at  hifeh  sliding  speeds  has  been  reported  In  only  a  few 
references,  preferably  with  metal  to  metal  contacts,  and  in  each  case 
some  evidence  or  theory  has  been  related  to  a  melt  lubrication.  The 
subject  was  stimulated  by  World  War  II  in  an  attempt  to  understand 
projectile  rotating  band  friction.  Herzfeld  and  Kbsson  [41],  in  a 
"Confidential"  1953  report  recently  unclassified,  cited  military 
evidence  of  metal  band  melt  for  projectiles  reaching  speeds  up  to  6600 
ft/sec.  They  proposed  a  hydrodynamic  slider  bearing  model  of  the  band 
riding  on  its  own  melt.  The  melt  viscosity  was  assumed  to  be  constant, 
and  a  simple  lumped  energy  balance  was  used  to' establish  the  melt 
thickness.  A  feature  of  their  model  was  the  heat  loss  tc  the  barrel  by 
transient  heat  conduction.  Film  thicknesses  on  the  order  of  10“5  to 
10“7  cm  were  predicted. 

Wilson  [42],  unaware  of  the  above  work  which  was  not  published  in 
the  open  literature,  recently  proposed  a  similar  one-dimensional 
hydrodynamic  sliding  model  for  melt  lubrication.  Again,  the  fluid 
temperature  was  assumed  to  be  constant  at  the  melt  temperature  and  no 
attempt  was' made  to  account  for  heat  transfer  from  the  melt  region. 
However,  Wilson  did  point  out  that  the  results  differ  widely  depending 
on  whether  the  slider  melts  or  its  stationary  track  melts.  Figure  14.' 
Bicego,  et  al.  [43]  modified  Wilson's  work  to  account  for  heat 
conduction  through  the  melt  film  due  to  a  track  teaqierature  higher  than 


the  melting  slider.  In  ,the  case  of  an  ice  slider,  a  small  temperature 
difference  (.2°C)  was  sufficient  to  make  conduction  prevail  over 
viscous  heating.  Results  were  obtained  in  terms  of  the  track-slider 
temperature  difference,  thus  avoiding  the  actual  heat  conduction 
problem. 

One  of  the  first  published  reports  on  high  speed  friction  was 
carried  out  by  Johnson,  et  al.  [44].  A  steel  ball  was  slid  on  a  steel 
disk  up  to  speeds  of  110  ft/sec.-  The  coefiicient  of  friction  decreased 
from  0.55  to  0.25  for  increasing  sliding  speed.  Amontons*  law  (u 
independent  of  load)  was  obeyed  up  to  and  including  load3  with  contact 
stresses  of  255,000  psi. 

Perhaps  the  most  comprehensive  experimental  study  of  metal 
friction  and  wear  at  high  sliding  speeds  was  done  at  the  Franklin 
Institute  from  1946  to  1956  by  Clark,  Morsell,  and  Shugarts.  There  was 
no  publication  of  this  work , in  the  open  literature  because  it  was 
classified  during  that  time.  Montgomery  [9]  has  recently  collected  all 
of  this  data.  A  pin-disk  machine  was  used  to  slide  various  metal 
specimens  against  a  steel  disk.  The  coefficient  of  friction  decreased 
with  increasing  values  of  pressure  x1  velocity  for  all  tested  metals.  A 
typical  low  coefficient  of  friction  was  approximately  0.2.  One 
reported  run  with  a  nylon  pin  gave  y  -  0.10  at  900  ,ft/3ec  and  4100  psi. 
Montgomery  found  that  the  pin  wear  rate  correlated  directly  with  the 
reciprocal  of  the  material  absolute  melting  point  (Figure  15). 

Montgomery  [10][11j  also  reported  on  experimental  work  to 
determine  the  friction  coefficient  of  projectile  rotating  bands  made 
from  gilding  metal  (90  CU;  10  Z^).  Measurements  of  gas  propellant 
pressure,  projectile  acceleration,  and  the  band  normal  contact  pressure 
on  rounds  fired  in  a  155  mm  howitzer  were  used  to  calculate  the  . 
friction  coefficient.  The  results  were  compared  with  similar  data  from 
the  pin-disk  experiments  of  the  Franklin  Institute  (Figure  16).  The 
friction  coefficient  for  rotating  bands  drops  quickly  to  a  steady-state 
value  of  0.02  at  a  pressure  x  velocity  of  approximately  0.8  x  10^ 
(psiKfps).  This  value  is  an  order  of  magnitude  below  the  pin-disk 
value.  Montgomery  [11]  attributed  the  difference  to  the  3ize  effect 
inherent  in  the  hydrodynamic  sliding  model  of  melt  as  proposed  by 
Wilson. 


Figure  16.  Coefficient  of  Friction  as  a  Function  of  the  Product  of  Bearing  Pre 
Sliding  Velocity  for  a  Projectile  with  Gilding  Metal  Rotating  Bands 


In  early  experiments  using  the  Franklin  Institute  pin-disk 
machine,  Sternlicht  and  Apkarian  [45]  measured  the  friction  and  wear  of 
several  metal  pins  (Mo,  W,  Cr,  Cu,  Ag,  Al)  on  a  steel  disk  with  speeds 
up  to  2000  ft/sec.  Again,  the  coefficient  of  friction  decreased  with 
an  increase  In  velocity.  They  measured  the  electrical  contact 
resistance  across  the  interface  of  the  rubbing  metals.'  An  increase  in 
contact  resistance  with  velocity  was  attributed  to  the  developing  ' 
molten  liquid  layer  at  the  interface.  In  an  application  of  Reynolds' 
equation  to  melt  lubrication,  iternlicht  and  Apkarian  proposed  the 
thermal  wedge  [461  (increasing  film  temperature  gradient  in  the 
direction  of  motion]  as  the  mechanism  to  support  a  normal  force  between 
parallel  surfaces  moving  with  relative  velocity.  Based  on  the  solution 
of  this  problem,  the  investigators  found  that  the  calculated  pin  wear, 
due  to  end  leakage  of  the  molten  film,  differed  from  the  measured  wear 
within  an  order  of  100J.  Neither  the  basic  equations  nor  the  solution 
details  were  offered,  in  the  paper. 

Bowden  and  Freitag  [47]  studied  high  speed  friction  by  spinning  a 
steel  ball  in  a  magnetic  field  and  measuring  the  deceleration  when  it 
was  brought  in  contact  with  a  metal  or  diamond  surface.  Melting 
occurred  at  the  region  of  contact.  An  analysis  of  the  heat  flow  into 
the  specimen  showed  that  the  area  of  intimate  contact  was  very  small 
compared  to  the  apparent  area.  . 

Bowden  and  Persson  [48]  contihued  the  spinning  ball  experiments  of 
Bowden  and  Freitas  and  presented  ail  in-depth  study  of  melt  with  high 
speed  friction.'  Figure  17  shows  the  results  with  a’  5-kg  load  for  a 

m.p.  65°C),  bismuth  (m.p.  271°C), 

),  copper  (m.p.  1080°C)  and  3teel 
s  of  friction  shifted  to  higher 


steel  ball  sliding  on  Woods  alloy 
tin  (m.p.  232°C)  ,  lead  (m.p.  328JC 
(m.p.  1500°C).  Minimum  coefficient 
sliding  velocities  as  the  metal  mel 
cases  the  minimum  coefficient  of  fr 


ting  point  increased.  In  several 
iction  was  well  below  a  value  of 


0.1.  A  substantial  increase  in  wear  occurred  at  sliding  velocities 
corresponding  to  the  minimum  p.  Examination  of  worn  surfaces  revealed 


large  formations  of  splashed-out  so 
edge  as  well  as  small  globules  of  S' 


lidified  material  near  the  leading 
olidified  melt  on  the  wear  track. 


This  evidence  was  much  less  pronounced  on  the  high  melting  copper  and 
steel.  Using  the  work  of  Landali  [4<i],  who  determined  the  propagation 


of  a  melt-solid  interface,  Bowden  and  Persson  crudely  calculated  the 
time  for  a  steady  state  melt  condition  to  be  reached,  approximately 
10"5  sec.  This  number  was  much  smaller  than  the  ball  contact  time. 
Figure  18  shows  the  results  for  several  non-metals:  silver  nitrate 
(m.p.  ?10°C),  nylon  (m.p.  265°C),  teflon  (decomposes  at  *iOO°C), 
butadiene  styrene  copolymer  rubber,  and  glass  (m.p.  800°C).  Both 
silver  nitrate  and  glass  developed  the  typical  wall  of  molten  material 
at  the  leaving  edge.  Teflon ’ displayed  a  rough  hide-lika  texture  at 
high  speeds  denoting  uneven  decomposition.  Molten  nylon  was  splashed 
out  of  the  wear  mark;  however,  the  surface  had  a  large  number  of  pits 
or  dimples,  and  the  extruded  material  contained  bubbles,  suggesting 
gaseous  matter  was  produced.  The  rubber  also  had  a  foamy  appearance  to 
the  molten  debris.  Another  approximate  calculation  by  Bowden  and 
Persson  indicated  that  steep  temperature  gradients  existed  across  the 

molten  nylon  film  with  film  temperatures  exceeding  800°C.  Metal  film 

temperatures  probably  were  only  a  few  degrees  above  the  melting 
temperature  even  at  the  higher  speeds.  This  difference  was  attributed 
to  large  melt  viscsosities  and  low  thermal  conductivities  of  polymers. 

Miller  [50]  investigated  the  surface  of  various  carbides  and 
oxides  with  high  melting  points  after  they  were  subjected  to  the  same 
spinning  steel  ball  apparatus  as  above.  Wear  from  a  network  of  surface 
creeks  correlated  with  their  lack  of  resistance  to  thermal  shock. 

Carignan  and  Robinowicz  [51]  have  carried  out  sliding  tests  of 

soft  metal  and  nonmetal  pins  on  a  rotating  steel  disk  at  speeds1  up  to 
150  m/sec.  The  friction  data  was  in  approximate  agreement  with  those 
obtained  by  others  using  the  pin-disk  geometry. ' 

7.  DISCUSSION 

A  review  of  the  literature  indicates  the  following  scenerio  for 
the  effect  of  sliding  velocity  on  the  coefficient  of  friction.  ■  For  low 
sliding  velocities  the  coefficient  of  friction  increases  with  speed, 
presumably  caused  by  an  increased  contact  area.  At  some  critical 
velocity,  on  the  order  of  several  feet  per  second  for  polymers,  the 
friction  coefficient  peaks  and  begins  to  decay  with  increasing  speed. 
Evidence  Is  overwhelming  that  material  melting  occurs  at  the  asperity 
.tips  for  these  relatively  low  sliding  velocities.  Melting  increases 


with  sliding  speed  until  the  entire  apparent  contact  area  is  a  layer  of 
melted  material.  The  speed  at  which  the  full  melt  layer  exists  is  not 
established,  but  it  does  depend  on  the  load.  The  data  from  several 
sources  suggest  that  several  hundred  feet  per  second  is  sufficient  for 
a  full  melt  of  polymer  material  at  light  loads  (pressures  much  less 
than  the  yield  stress).  Although  several  authors  in  the  early  70’3 
discussed  the  importance  of  apparent  pressure  as  a  parameter  in 
friction  work  with  polymers,  load  continues  to  be  used  instead.  Little 
work  was  done  at  high  loads  with  the  exception  of  internal  ballistics. 
Montgomery  found  that  silding  metal  produced  a  full  melt  layer  at  PV  - 
5  x  10^  (psi)Ifps).  For  measured  contact  stresses  of  50,000  psi,  this 
limit  giyes  a  sliding  velocity  of  100  ft/sec.  Polymers  have  much  lower 
thermal  conductivities  and  melting  points  than  metals. 

Two  recent  references  [42][43]  focused  on  a  melt  lubrication 

t 

model  to  explain  the  low  friction  of  ice.  This  model  relied  on  the 
hydrodynamic  3lider  bearing  concept.  Montgomery  noted  that  the 
coefficient  of  friction  for  gilding  metal  from  internal  ballistics  (m  - 
0.02)'  tests  was  milch  lower  than  values  found  for  pin-di3k  experiments 
(u  >  0.1).  It  is  significant  that  most  experimental  work  on  friction 
was  accomplished  with  pin-disk  machines,  all  pin  diameters  being  nearly 
1/10  inch.  All  reported  coefficients  of  friction  on  these  machines 
were  greater  than  0.1  at  high  speeds.  The  one  exception  was  the  work  , 
done  by  Bowden  and  Persson  [48]  on  spinning  steel  balls.  Here, 
friction  coefficients  well  below  0.1  were  recorded. 

One  difficulty  with  the  slider  bearing  concept  is  the  need  for  a 
lubricant  at  the  front  edge  of  the  slider.  Wilson  [42]  was  well  aware 
of  this  problem  and  he  recommerded  a  minimum  load  before  melt 
lubrication  can  occur.  Montgomery  seized  upon  the  length  factor  of 
this  minimum  load  to  explain  the  difference  in  friction  coefficients 
between  sliding  projectile  bands  and  pin-disk  experiments  at  comparable 
conditions.  This  writer  will  propose  an  alternate  solution  to  the  melt 
lubrication  problem  below  which  circumvents  the  lubricant  problem; 
however,  there  is  a  better  explanation  for  the  high  friction 
coefficients  achieyed  with  pin-disk  experiments.  If  the  friction 
coefficient  is  calculated  from  3lider  bearing  theory  [46],  it's  value 


is  proportional  to  a  functional  of  the  length  to  width  ratio  which 
characterizes  rotating  bands,  when  compared  to  the  ratio  of  one  for 
pins. 
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SECTION  III 


LITERATURE  REVIEW  OF  RHEOLOGY 

1.  NEWTONIAN  BEHAVIOR  (OILS) 

For  Newtonian  fluids  the  shear  stress  is  linearly  related  to  the 
strain  rate,  i.e., 


The  absolute  viscosity  u  is  generally  dependent  on  pressure  and 
temperature.  Early  measurements  made  at  Harvard  and  sponsored  by  ASME 
[52]  show  this  dependency  over  a  large  range.  They  suggested  an 
empirical  relationship 

ln(— )  -  c  e®/T  {!  ♦  ap  -  fl[ap  -  e~“P)]}  (18) 

w° 

The  constants  y0,  c,  B,  o,  8  vary  for  each  oil. 

At  atmospheric  pressure  the  variation  of  viscosity  with 
temperature  is  given  by 

ln(— )  -  c  e6/T  (,9) 

Ho 

The  relationship  between  temperature  and  either  Saybolt  seconds  or 
kinematic  viscosity  may  be  plotted  aa  a  straight  line  on  charts  with 
special  scales,  known  as  the  A STM  Standard  Viscosity  Temperature  Charts 
for  Liquid  Petroleum  Products  D341.  However,  straight  lines  for  . 
absolute  (dynamic)  viscosity  may  be  plotted  with  minor  error  as  shown 
in  Figure  19. 

In  most  hydrodynamic  theory  the  effect  of  pressure  is  not  t;<ken 
into  account.  This  effect  can  be  ignored  up  to  a  fluid  pressure  of 
1000  psi,  but  at  pressures  of  5000  psi  ordinary  lubricating  oils  double 
their  viscosity.  In  1&93  Barus  established  an  empirical  equation  to 
describe  isothermal  viscosity-pressure  relationship  for  a  given 
liquid: 

—  -  e°P  (20) 

Ho 
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Figure  19.  Dynamic  or  Absolute  Viscosity  vs.  Temperature  for 
SAE-Nuabered  Oil 


Dynamic  vifeoiity  o,  fl»P»  *  I  (IctmipoiK)  Of  10  *  MPa  ■  « 


Eyring  and  co-workers  [53]  have  derived  a  similar  relationship  based 
upon  molecular  free  volume  mobility  theory.  Generally  1/a  is 
approximately  equal  to  7000  psi  for  most  mineral  oils,  resins,  and 
polymer  blends.  An  early  summary  on  pressure-viscosity  work  was 
published  by.  ASME  [54],  see  Figure  20. 

Interest  in  viscosity-pressure  correlations  has  paralleled 
interest  in  elastohydrodynamics  (EHD).  .Equation  (20)  above 
overestimates  the  viscosity  at  very  high  pressures.  Cameron  [55]  has 
proposed 

—  -  (1  +  Cp)n  (21 ) 

Ho 

Other  correlations  are  discussed  by  So  and  Klaus  [56]  and  Johnson  [57]. 
Allen,  C.  W.,  et  al.  [58]  reported  that  spinning  ball  experiments 
showed  poor  agreement  between  load  and'  torque  using  equation  (20). 

They  proposed  a  modified  exponential  form  (Figure  21): 


—  -  eaP  PiPg, 

Ho 

-  e<*Pm  ♦  0(p  -  pm)  p  >  pm 


(22) 


A  similar  model  was  used  successfully  by  Cheng  [59].  However,  Johnson 
[60]  argues  that  the  reduced  viscsosity  at  very  high  pressures  is 
caused  by  the  shear  stress  limiting  feature  of  the  oil  at  these 
pressures.  Mon-Newtonian  fluid  behavior  is  now  discussed. 


2.  NON-NEWTONIAN  3EHAVI0R  (OILS)' 

Non-Newtonian  fluid  behavior  is  characterized  by  a  viscosity  which 
varies  with  rate  of  3hear.  Polymer-th tokened  oils  (oils  containing  VI 
improvers)  usually  suffer  viscosity  losses  when  subjected  to  high  shear 
rates.  These  losses  can  be  either  permanent  (due  to  a  mechanical  or 
chemical  breakdown  of  the  larger  molecules)  over  a  period  of  time,  or 
it  can  be  temporary,  recovering  the  original  viscosity  when  the  stress 
is  removed.  We  are  concerned  with  the  latter. 

The  American  Petroleum  Institute  (API)  in  conjunction  with  the 
ASTM  conducted  an  investigation  on  viscsosity  variation  with  shear 
rates,  reported  in  ASTM  Special  Technical  Publication  No.  Ill: 
"Symposium  on  Methods  of  Measuring  Viscosity  at  High  Rates  of  Shear." 
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Figure  20.  Viscosity-Pressure  Relationships  for 

Paraffinic  and  Naphthenic  Oils  at  Three 
Temperatures 
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Figure  21.  Theoretical  Pressure-Viscosity  Relation  for 
Synthetic  Paraffinic  Oil  at  83°F 


Two  polymer  blended  oils  underwent  a  30-40  percent  temporary  viscosity 
loss  at  shear  rates  of  50,000  to  1,000,000  sec"1 —comparable  to  rates 
found  in  automobile  engines.  Klaus  and  Fenske  [61]  reported  on 
viscosity-shear  characteristics  of  lubricants  and  related  the  results 
to  hydraulic  systems.  West  and  Selby  [62]  studied  the  effect  pf 
multigraded  polymer  blended  engine  oils  and  found  permanent  viscosity 
losses.  In  January  of  1964,  a  symposium  was  held  by  ASTM  and  published 
in  ASTM  Special  Technical  Publication  No.  382:  "The  Effects  of  Polymer 
Degradation  on  Flow  Properties  of  Fluids  and  Lubricants  Containing 
Polymers." 

In  the  mathematical  description  of  non-Newtonian  fluids  it  is 
necessary  to  retain  the  basic  momentum .equation  in  lubrication: 


3p  jk 
3x  3y 


(23) 


Flow  properties  of  the  lubricant  are  expressed  in  one  of  two  general 
forms : 


3u 

3y 


F-|  (  t  ) 


(24) 


or  the  inverse 

t  »  F2(3u/3y)  •  (25) 

Greases  have  been  modeled  as  Bingham  plastics  by  Milne  [63H64]: 


t  »  ±tq  + 


(26) 


in  which  a  yield  value  must  be  reached  before  flow  takes  place. 

Examples  can  be'  founo  in  Pinkus  and  Sternl icht  T65] .  If  a  liquid  is 
stressed  rapidly  enough,  it  will  display  an  elastic  response  to  stress; 
thus,  the  duration  of  the  stress  becomes  a  parameter  irt  describing  the 
properties  of  such  fluids.  Milne  [66]  has  applied  a  Maxwell  model  of  a 
lubricant  (linear  viscoelastic  fluid)  to  the  theoretical  treatment  of 
slider  bearings: 


3u  1  dt  t_ 
3y  G  dt  p 


(27) 


were  G  is  the  elastic  modulus  and  the  viscous  dissipation  term  (second 
on  right)  is  Newtonian.  If  conditions  are  steady,  the  time  derivative 
may  be  replaced  by  the  space  derivative  and 

.1  •  (»> 

3y  G  9x  y 

Non-Newtonian  effects  can. also  be  expressed  in  terms  of  non-linear 
viscous  dissipation  and/or  shear  limiting  effects.  Smith  [67]  was  the 
first  to  suggest  that  lubricants  under  high  pressure  contacts  display 
characteristics  that  are  solid-like.  Johnson  and  Cameron  [68]  carried 
out  experiments  using  cylindrical  disks  in  combined  rolling  and  sliding 
for  high  contact  pressures.  They  reported  curves  similar  v-o  Figure  21; 
however,  Johnson  [60]  speculated  that  the  lubricant  behaved  a3  a 
granular  solid,  displaying  a  critical  shear  stress  which  can  not  be 
exceeded,  rather  than  possessing  a  reduced  viscosity  at  uhe  high 
pressures. 

a.  Glass  Transition 

When  lubricants  are  cooled,  they  do  not  solidify  into  a 
crystalline  solid  at  some  fixed  temperature.  Instead  their  viscosity 
steadily  increases  until  the  material  exhibits  a  solid-like  behavior 
referred  to  as  a  glass  (amorphous  ductile  solid).  Although  the 
process  is  continuous,  properties  such  as  specific  volume  or  speed  of 
sound  show  a  marked  change  at  the  glass  transition  temperature.  The 
glass  transition  temperature  is  not  fixed.  Since  the  transition  is  a 
viscoelastic  phenomenon,  it  responds  in  a  glassy  (elastic)  way  at 
ambient  pressures  and  temperatures  a?  in  the  early  ASTM  work  above.  A 
major  method  to  test  the  viscoelasticity  of  oils  to  strain  rates  is  the 
use  of  oscillatory  shear  at  frequencies  to  78  MHz.  The  now  classic 
paper  in  this  field  is  that  of  Barlow  and  Lamb  1691.  Major  ' 
contributions  have  been  made  by  Dyson  [70][7.1]  to  incorporate  the  above 
work  into  a  continuous  shear  mechanism  of  EHD.  Differences  between 
oscillatory  shear  and  EHD  experiments  to  extract  the  elastic  shear 
modulus  is  attributed  to  relaxation  times  by  Johnson  [72].  A  good 
summary  of  the  oscillatory  shear  approach  to  rheology  is  given  by 
Hutton  [73]. 
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Hirst  and  Moore  [74]  have  shown  that  a  viscosity  in  the  EHD 
contact  region  of  1C>5  pa-s  (14.5  lb-sec/in^)  is  a  good  working 
criterion  for  the  onset  of  an  elastic  response.  If  the  glass 
transition  is  a  consequence  of  the  viscosity  reaching  a  critical  value, 
the  transition  could  equally  as  well  be  brought  about  by  inoreasing  the  r 

pressure  isothermally  as  by  cooling  at  constant  pressure.  Johnson  and 
Roberts  [75]  revealed  viscoelastic  behavior  in  a  point  contact  disk 
machine  at  contact  pressures  greater  than  0.5'  Gpa  (72,500  psi). 

Sliding  speed  (strain,  rate)  was  very  low.  This  aspect  of  glass 
transition  has  been  studied  by  Alsaad.et  al.  [76].  They  observed  that 
glassy  states  are  very  likely  for  most  lubricants  in  EHD  contacts  and 
occur  at  lower  pressures  for  sliding  as  opposed  to  rolling  contacts. 

Also,  the  glass  transition  temperature  rises  witn  increasing  pressure, 

•  Figure  22. 


b.  Johnson  Rheology  Model 

The  impetus  for  viscoelastic  studies  of  lubricants  in  EHD 
contacts  was  the  tracticr  problem  using  disk  machines  (two  parallel 
rollers  with  different  peripheral  speeds.  Typical  traction  curves  [77] 
are  shown  in  Figure  23.  In  general  there  are  three  regions:  (i)  a 
linear  Newtonian  region,  (ii)  a  non-linear  region  where  the  traction 
approaches  a  constant  value  independent  of  the  rate  of  shear  (slip), 
and  (iii)  a  thermal  region.  The  traction  is  dominated  by  the 
non-linear  region.  Hirst  and  Moore  [78]  first  presented  a  non-linear 
viscous  flow  model  for  high  pressure  EHD  based  upon  the  Eyring  thermal 
activation  theory.  Expanding  on  this  idea  Johnson  and  Tevaarwerk  [79] 
developed  a  non-linear  Maxwell  model  described  by 


du  1  dt 
dy  .  G  dt 


F(-t) 


(29) 


where 


F(x)  =  r2’  sinh(x/T0)  •  (30) 

u  ' , 

and  t0  is  an  experimental  constant.  At  small  shear  stresses  F(t) 
reduces  to  ,t/y — a  Newtonian  fluid.  It  is  important  to  realize  that  G 
and  r0  depend  upon  temperature  and  pressure  as  does  the  viscosity.  The 
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Figure  22.  Heuristic  Estimates  of  the  Relationship  Between 
Conditions  in  an  EHD  Contact  and  Glass-Liquid  , 
Phase  Diagrams  of  Some  Lubricants 
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Figure  23.  Typical  Traction  Curves  Measured  on  a 
Two-Disk  Machine  in  Line  Contact  at 
Varying  Mean  Contact  Pressures 


success  of  this  model  is  shown  in  Figure  24  which  includes  data  from  a 
disk  machine.  The  dimensionless  number  D  CDeborah  number)  «  yU/Gb  ■ 
where  b  is  the  contact  radius.' 

Eyring's  equation  for  non-linear  shear  rates  may  be  written 


du 

dy 


A  exp - ! — -  sinh 

L  kT 


lill 
kT  J 


(31 ) 


where  A  is  a  constant,  E  is  the  activation  energy  at  atmospheric 
pressure,  T  is  the  absolute  temperature,  k  Is  the  Boltzmann's  constant, 
vi  and  V2  are  activation  volumes  for  pressure  and  shear,  respectively . 
When  v2t  >>  kT 


du 

dy 


1 

-  -  A  exp 


v2 i-vip-E 
kT 


(32) 


or  comparing  with  equation  (30) 
x0  -  kT/v2 


In  the  high  strain  rate  region  equation  (32)  implies  that,  at  constant 
strain  rate. 


♦  constant 


(33) 


Such  a  relationship  has  been  recently  found  for  the  shear  of  solid, 
amorphous  polymers  [80][8l ][82][83][84]. 

,c.  Winer  Rheology  Model 

_  Experiments  have  been  undertaken  by  , Bair  and  Winer  [8.5]  to 

measure  the  shear  stress  versus  shear  rates  (low  rates  10-1*  s-1  to  high 
rates  10.2  s-1)  for  various  lubricants  under  high  pressure  contacts,  see 
Figure  25  for  example,  the  data  lies  primarily  in  the  high  shear 
region  which  exemplifies  the  limiting  shear  stress.  This  data  is 
compatible  with  Johnson  above  who  acknowledged  a  limiting  shear  stress 
by  advocating  a  Prandtl-Reuss  model,  (elastic-plastic  solid)  at  very 
high  rates  of  strain.  Bair  and  Winer  were  able  to  collapse  the  data 
for  several  lubricants  by  means  of  the  non-linear  viscous  model 


F(t)  =  -  — ^  Hn(l  -  — )  .  '  •  (3*0 

u  tl 

where  t ^  is  the  limiting  shear  stress.  Again  depends  on  pressure 
and  temperature.  Recent  experimental  data  by  Bair  and  Winer  [36]  for 
twelve  lubripants  give  pressure  and  temperature  effects  on  viscosity, 
elastic  modulus,  and  the  limiting  shear  stress.  See  also  Houpert.et 
al.  [87]. 

3.  POLYMERS 

Oils  under  the  extreme  pressures  of  elastohydrodynamie  (EHD) 
contacts  thicken  to  the  point  of  becoming  plastic-like.  Thus,  their 
behavior  is  very  similar  to  polymer  melts.  The  viscosity  of  polymers 
have  been  described  by  expirical  relations  for  over  a  half  of  a  century 
while  the  description  of  oils  under  EHD  contacts  is  very  recent. 

Polymer  rheology  is  now  reviewed. 

Eyrings  theory  [88],  outlined  for  oils  in  Section  2b.,  was 
proposed  to  explain  the  shear  behavior  of  polymers  as  early  as  1936. 

Let  the  shear  rate  du/dy  =  Y.  Then  rearranging  equation  (3T) 


T 

U  -  ~ 


A  expt^  P/T)  +  (a2/T)]  sinhiajr/T) 


If  temperature  changes  are  not  too  large  when  pressure  and  shear  stress 
changes  occur, 

v  -  fi(M)f2(T)f3(P),  fij(x)  (36) 

where  M  is  the  molecular  weight.  It  has  been  found  experimentally  that 
most  unbranched  linear  polymers  having  a  molecular  weight  in  excess  of 
the  critical  Mc  obey  the  above  rule  of  logarithmic  additivity  of 
viscosity.  A  linear  polymer  is  a  polymer  in  which  the  molecules  are 
linked  together  in  the  form  of  chains  with  little  side  links  or 
branches  such  as  high  density  polyethylene. 

At ' low  shear  stress  (or  shear'  rates)  equation  (35)  becomes 
Newtonian  but  with  increasing  shear  stress  the  polymer  melt  viscosity 
decreases  in  a  characteristic  way,  giving  rise,  to  what  is  known  as  a 
non-Newtonian  flow.  A  number  of  empirical  equations  have  been  proposed 
to' describe  this  behavior.  The  best  known  is  the  power  law  model  of 
Ostwald-de  Waele  (1923): 


--3 


Figure  26  illustrates  this  law  very  well  [89].  Thus, 


u  -  cyN-1  (38) 

or 

„  .  c.t(M-1)/H  (39) 

where  M  is  called  the  "flow  index"  of  the  polymer.  From  equations  (38) 
and  (39)  it  is  not  surprising  that  the  shear  stress  functional  is  often 
replaced  by  a  shear  rate  functional  in  equation  (36)  for  linear 
polymers. 

Of  course  the  power  law  can  not  adequately  represent  the 
ncn-Newtonian  flow  behavior  over  a  wide  shear  rate  range  nor  can  it 
predict  the  transition . to  the  low-shear  limiting  Newtonian  viscosity. 
Mendelson  [90]  proposed  the  following  for  a  good, fit  to  experimental 
data  at  low  shear  rates: 


log  t  -  CQ  ♦  Ci  log  Y  ♦  C2(log  Y)2  (40) 

Parrini  et  al.  [91]  investigated  the  melt  rheology  of  nylon  -  6,'  nylon 
-  6,  10,  and  nylon  -  11  by  different  techniques  and  arrived  at  the 
following  relationship: 


log  T 


-2.88  ♦  2.18  10g(u0Y)  -  0.121[log(uoY)]2 


(41) 


which  is  valid  for  shear  stresses  upto  1500  psi.  A  recent  correlation 
by  Mount  and  Chung  [92],  which  includes  temperature,  is  shown  in  Table 
1  for  several  polymers. 

Recognizing  the  limiting  viscosity  at  high  shear  rates  Cross  [93] 
proposed  a'  semlemplrlcal  relationship: 


♦  -Mo.  :-.»»»■ . 

1  ♦  (AY)“ 


where  A  is  a  temperature  dependent  constant  and 


(42) 


m  •  (M|f/Mg)0.2 

However,  the  relationship  is  a  four-parameter  (unknown)  equation. 
Other  empirical  relations  include 


TABLE  1.  PHYSICAL  PROPERTY  DATA  OF  POLYMERS  USED  IN  COMPUTER  S 


ft 


vy 

8 

£ 

a 


Bueche  and  Harding  C94] 


»  [1  *  O.6(XY)0.75]-1 
Wo 


where  X  -  12p0M/it2pkT 
Vinogradov  [95] 


U 

Wo  L 


.  0.355  .  0.71  |  -1 

-  1 1  •  o.«(i^)  •  (JSl) 
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(44) 

(45) 


where  C  »  1.9  psi 
Graeasley  [96] 


wo  L 


♦  1.916(t0V/2) 


>.aoa*|  -i 


(46) 


The  above  empirical  fits  do  not  include  the  high  shear  rate  asyptotes. 
However,  considering  the  denominator  powers  of  the  shear  rate  are 
approaching  one, 

t®  -  Praia,  »  po/(0.6  X)  eq.  (43) 

*  u0/t0  eq.  (46) 

i.e.,  the  shear  stress  approaches  a  limit  at  high  shear  rates  as  with 
lubricating  oils. 

The  temperature  function  in  equation  (36)  is  generally  givein  by 
the  Eyring  relation  for  both  lubricating  oil3  and  polymers: 

f 2 ( T )  -  Aeu/kT  (47)  ' 

or 

,1  1 

*  *0 


Wo 

Cog3well  [97]  lists  the  relative  fluidity  index 

'  RFI  .  iLilllc - 

y  at  T  +  10°C 

for  a  number  of  polymers,  .see  Table  2: 


(48) 


(49) 
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It  is  straight  forward  to  show  that 
in(RFI) 

8  ”  - - : - 

i  1 

T  T  1 0°K 


(50) 


The  effect  of  increased  pressure  on  viscosity  is  similar  to  a 
decrease  in  temperature  i.e.,  reducing  free  volume  and  molecular 
mobility  so  leading  to  an  increase  in  viscosity.  The  study  of  the 
influence  of  pressure  on  viscosity  has  attracted  few  experimentalists 
or  detailed  reviews.  This  is  because  the  primary  interest  in  polymer 
properties  is  extrusion  processes  where  the  pressures  of  several 
hundred  psi  are  insufficient  to  influence  viscosity.  A3  with 
lubricating  oils 


u 

Ho 


ap  + 
»  e 


Since 


^  ap 

3P  Uo 


(— )  dT 
dT  u0 


(51 ) 


or 

(~)  -  oTq2/8  (52) 

dp  u0 

Cogswell  [97]  has  given  a  table  of  the  pressure/temperature  equivalence 
at  constant  viscosity  as  well  as  constant  entropy  to  suggest  a  guide  if 
only  the  thermodynamic  data  is  available.  The  table  is  reproduced 
below.  With  Table  3  and  equation  (52),,  the  pressure  coefficient  can  be 
found  for  a  number  of  polymers.  , 


4.  ADOPTED  MODEL 

Both  the  lubricating  oils  (under  high  pressures)  and  the  polymers 
have  similar  properties  and  have  similar  mathematical  descriptions. 

They  show  elastic  behavior  (memory)  during  short  time  scales  so  that 
equation  (29)  is  the  fundamental  starting  point  for  a  model.  A  classic 
text  on  the  viscoelastic  properties  of  melt  polymers  is  written  by 
Ferry  [98].  He  gives  an  expression  which  relates  the  elastic  shear 
modulus  G  to  the  molecular  weight  distribution  of  the  polymer:  • 


TABLE  2..  RELATIVE  FLUIDITY  INDEX  (RFI) 
FOR  AN  INCREASE  IN  TEMPERATURE  OF  10 °C 


Polymer 

Test 

temperature 

ra 

RFI* 

Branched  polvethvlcne 

150 

1-35 

2IX) 

13 

250 

1-25 

Linear  polyethylene 

200 

12 

Polypropylene 

200 

1-2 

Polystyrene 

2(KI 

1-7 

Polvmethvl  methacrvlate 

21X1 

2-5 

' 

250  ' 

10 

Polycarbonate 

250 

15 

Polyethersulphone 

350 

1-5 

Nylon  6:6 

275 

1-35 

Polyethvle"  ••  terephthslate 

£75 

F-35 

Polyacef 

21X1 

1-2 

Polyv  aylchloride 

21X1 

3-0+ 

•  RFI  =  — viu-osh\  .ii  7  f  „  cons|;lni  MIVv.. 

viscosity  ut(7  +  10)  C  .  , 

t  In  the  case  of  pulyvinjlchloride  the  meft’  is  known  to  he  slicbtly 


crystalline. 


TABLE  3.  RATIO  OF  VISCOSITY  TO  ENTROPY 
OF  VARIOUS  POLYMERS 


Polymer 

(ATIAP). 

(•C/Nm-=) 

(dTIdPU 

(•C/Nm-i) 

Ratio.  17/S 

Low-density  polyethylene 

5-3  x  lO-’ 

1-6  x  I0-7 

3-3 

High-density  polyethylene 

4-2 

1-5 

2-8 

Polypropylene 

8-6 

2-2 

4-0 

Polystyrene 

4-0 

1-5 

2-7 

Poly( methyl  methacrylate) 

3-3 

1-2 

2-8  ' 

PVC 

3-1 

11 

2-8 

Nylon 

3-2 

1-2 

2-7 

Silicone 

6-7 

19 

3-5 

Acetal  copolymer 

51 

,  1-4 

3-6 

Polyphenylene  oxide 

5-7 

— 

— 

Polycarbonate 

3-6 

— 

— 

Polyethersulphone 

6-7 

— 

— 

Averagi  val"e 

5-0  SD  1-7 

1-4  SD  0-4 

3-1  SD  <1*5 

\ 

\ 


1  0-4 
G  =  pRT 


(53) 


where  Mw  is  the  weight  average  molecular  weight  as  opposed  to  Fty,  the 
number  average,  and  R  is  the  gas  constant.  Generally,  the  modulus  is  ' 
on  the  order  of  10^  pa  (15  psi)  for  both  lubricating  oils  and  polymers 
Polymer  models  must  have  Newtonian  behavior  at  low  rate3  of  shear 
and  limiting  shear  stress  at  high  rates  of  shear.  Between  these  two 
limits  a  number  of  empirical  models  have  been  reviewed  above.  Often, 
there  is  little  to  choose  between  them  when  selecting  the  proper 
experimental  constant.  Since  rotating  bands  operate  in  very  high  shear 
rate  regions  where  limiting  shear  stresses  exist  near  the  moving 
surface,  it  is  proposed  to  use  the  more  direct  (and,  mathematically, 
the  least  complex)  low  and  high  shear  rate  asymptotes.  Thus, 


au 

dy 


1_(1) 

dt  G 


+  F(t) 


(5*) 


where 


or 


t  du 
’  tl  dy 


t  <  tl 

T  -  TL 
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SECTION  IV 


FRICTION  THEORY  FOR  METAL  ROTATING  BANDS 
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INTRODUCTION 


Evidence  of  material  melt  at  contact  points  of  rubbing  surfaces  is 
well  documented  by  Bowden  and  Tabor  [14].  That  full  melt  occurs  over 
the  surface  at  high  sliding  speeds  is  not  surprising.  Bowden  and 
Freitag  [47]  found  a  melt  over  the  contact  surface  of  metals  touched  by 
spinning  steel  balls.  Sternlicht  and  Apkarian  [45]  gave  evidence  of 
melt  with  the  Franklin  Institute  pin-disk  machine.  Montgomery  [9][10] 
proposed  a  full-film  melt  for  projectile  rotating  bands. 

Under  large  loads  full-film  melt  can  take  place  at  low  sliding 
speeds.  Tanaka  [40]  reported  a  thin  layer  of  melt  on  polymer  pins  when 
rubbed  with  steel  and  gla33  surfaces  at  speeds  to  2.5  m/sec.  Recently 
Clerico  [38]  found  the  presence  of  a  solidified  melt  on  rotating 
polyacetal  disks  at  relative  sliding  speeds  of  0.76  m/3ec. 

Theoretical  papers  on  melt  lubrication  are  limited.  Wilson  [42] 
proposed  a  hydrodynamic  slider  bearing  model  of  the  slider  riding  on 
its  own  melt.  The  fluid  temperature  wa3  assumed  to  be  constant  at  the 
melt  temperature,  and  no  consideration  was  given  to  heat  transfer  from 
the  melt  region.  Bicego,  et  al.  [43]  modified  Wilson's  work  to  include 
conduction  from  the  surface  (ice)  to  the  film.  The  results  were 
obtained  in  terms  of  the  track-slider  temperature  difference,  thus 
avoiding  the  actual  heat  conduction  problem.  There  is  no  a  priori 
reason  to  assume  that  non-parallel  surfaces  (hydrodynamic  effect)  are 
necessary  to  produce  load  support  for  melt  lubrication.  In  early 
experiments  using  the  Frarklin  Institute  pin-disk  machine,  Sternlicht 
and  Apkarian  [45]  measured  the  friction  and  wear  of  several  metal  pins 
against  a  high  speed  steel  disk.  In  an  application  of  Reynolds' 
equation  to  melt  lubrication,  they  proposed  the  thermal  wedge  as  the 
mechanism  to  support  a  normal  load  between  parallel  surfaces.  Neither 
the  basic  equations  nor  the  solution  details  were  offered  in  the  paper. 
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The  author  would  like  to  present  an  overlooked  mechanism  to 
provide  load  support  for  parallel  surfaces  undergoing  melt 
lubrication — the  mass  melt  itself  per  unit  area.  The  mechanism  is 
similar  to  the  load  support  achieved  with  porous  bearings. 

2.  ANALYSIS 

The  following  assumptions  are  now  made  but  verified  through  the 
course  of  the  section. 

1 )  The  fluid  is  laminar  and  incompressible. , 

2)  The  pressure  is  constant  across  the  film  and  average  values 
of  density  and  viscosity  are  employed. 

3)  Film  thickness  and  melt  penetration  into  the  solid  are  small 
compared  to  the  length  and  width  of  the  slider. 

;t;  Only  one  surface  is  melting.  Resolidification  occurs  after 
the  melt  flows  from  between  the  surfaces! 

5)  Quasi  steady-state  conditions  are  reached  In  the  film. 

The  basic  equations  for  thin  films  are  taken  from  Constant iqescu 


[10]. 


Momentum 

3P  3  ,  3u. 

(55) 

-  f  (,  £>  (56) 

3z  3y  3y 

Since  the  pressure  is  independent  of  y,  the  above  equations  may  be 
integrated  to  give  expressions  for  the  velocities  u,  w  between  the  two 
surfaces.  For  a  slider  moving  with  velocity  U  across  a  fixed  plane, 

u.l_i£(y2.yh).ffi  .  (57) 


u-k£  ',2  ■ s,n) 


where  h  is  the  film  thickness. 


Continuity 


3(puv  +  3 ( pw )  ^  3( pv) 
3x  3z  3y 


i  ?■  v  V 
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Whei.  the  continuity  equation  is  integrated  across  the  film,  the 
resu.tant  equation  is  called  Reynolds'  equation: 


1_  (S^l  i£\  +  i_  /  P^3  i£  \ 

3x  \  u  3*  /  3z  \  y  3z  / 


6U  —  (ph)  +  12[pv]h  (60)' 

ax  o 


Typically  in  such  derivations  the  term  pv  is  replaced  by  the  surface 
rate  of  approach,  p9h/3t,  and  the  film  is  referred  to  as  a  "squeeze 
film."  Previous  models  of  melt  lubrication  have  proposed  hydrodynamic 
action,  h(x),  or  thermal  wedge  action,  p(x),  to  account  for  film  load 
support.  However,  it  is  a  mistake  to  assume  that  the  last  term  in 
equation  (60)  is  zero  if  the  surfaces  do  not  .approach  each  other.  The 
term  should  be  interpreted  as  the  rate  of  mass  addition  per  unit  area 
from  the  melting  surface.  Then  for  parallel  surfaces  and  a  fluid  with 
constant  density  and  viscosity  equation  (60)  becomes 


V2p  -  12y  m/ph3 


(61) 


where  m  is  the  constant  mass  rate  of  melt  per  unit  area  for  the  slider 
or  fixed  surface.  A  solution,  analogous  to  the  squeeze  film  problem 
[98],  is  given  by 

48mt2y  •  (n-1)/2  fcosh(nirz/)l) 

P  "  Pa  "  ph3*3  nSl.3,5^15  "  |co3h(n*3/2) 

-  1  cos(nirx/t)^  (62) 

where  S  is  the  ratio  of  the  long  side  length  to  the  short  side  length 
l.  The  load  support  is  found  by  integrating  eauation  (62)  over  the 
ar~;  A. 

0  “  7  “  I" 1  *  ~5“  E,  ,  ,  n‘5tanh(nw8/2)l  (63) 

A  ph-J  l  *28  n-1,3,5  j 

or 

o  ■  ymt2{/ph3  (64) 

where  6  is  equal  to  one  for  large  8. 

The  energy  equation  is 


(65) 


kV2T  +  $ 


where  $  is  the  dissipation  function.  By  dimensional  reasoning  the  film 
is  quasi  steady-state  if  the  characteristic  time  of  the  process  tc  » 
h2/a.  Furthermore,  gradients  across  the  thin  film  are  much  greater 
than  gradients  in  the  film  plane,  and  the  latter  are  neglected  in  thin 
film  lubrication  [99].  Thus, 


(66) 


The  solution  is  found  by  substituting  equations  (57K58)  into  equatibn 
(66)  and  integrating  twice.  However,  a  simplification  can  be  made.  In 
all  pratical  cases  the  viscous  heating  due  to  the  squeeze;  velocity, 
which  arises  from  the  pressure  gradients,  is  negligible  when  compared 
to  the  heating  due  to  the  sliding  velocity.  Assuming 


h  3p 
2u  3x 


/(U/h) 


oh2 

ulu6 


(67) 


then 


T  '  Tf  *  (Ts  -  Tf)y/h  -  (y2  -  yh)  (68) 

where  T3 ,  Tf  are  the  surface  temperatures  of  the  slider  and  fixed 
surface  respectively. 

The  heat  flux  per  unit  area  out  of  the  film  can  be  established 
from  the  temperature  gradient  at  the  surface: 

qs  -  -k[dT/dy]h 

.  Qf  -  ♦k[dT/dy]0 

or 

9S  ’  k(Tf  -  Ts)/h  ♦  yU2/2h  -  (69) 

qf  -  -k(tf  -  Ts)/h  ♦  yU2/2h  (70) 

At  the  sc.m  time  the  heat  flux  into  the  surfaces  depends  upon 
whether  or  not  the  surface  is  melting.  Landau  [100]  has  solved  the 
transient  heat  conduction  problem  for  a  melting  solid  in  which  the  melt 
Is  continuously  removed.  He  showed' that  melting  cownences  at  time 
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tB  -  a[pHa/q]l/2  (71 ) 

where  a  is  the  diffusivity,  H  is  the  Latent  heat  of  fusion, 

a  =  ir1/2  c(Tm  -  T0)/2H  (72) 

and  Tm  is  the  melt  temperature.  For  rubbing  solids, 

•  q  «  foil  (rubbing  solids)  (73) 

where  a  is  the  load  per  unit  area,  and  f  is  the  friction  coefficient. 
For  highly  stressed  metal  surfaces  rubbing  at  moderate  speeds,  tm  is  on 
the  order  of  10-8  gec.  Landau  notes  that  a  steady-state  solution, 
i.e.,  a  state  characterized  by  the  migration  inward  at  a  constant 
velocity  of  a  fixed  temperature  distribution  in  the  solid,  exists  with 
the  heat  flux: 


qm  -  m[H  +  c(Tm  -  T0)]  (melt  surface)  (74) 


The  approximate  condition  at  which  this  steady-state  is  reached  is  t  * 
10  tm  for  ail.  With  smaller  values  (of  a)  the  time  is  longer.  At 
the  solid  surface  the  heat  flux  for  transient  conduction  in  a 
semi-infinite  solid  is 


q 


k’(T3.f  -  Tpj 
(ira't)'/2 


(75) 


Using  the  average  heat  flux  over  the  characteristic  contact  time  tc  - 
tc/U  where  the  characteristic  length  is  L  for  the  fixed  surface  and 
xQ  (the  sliding  distance)  for  the  slider, 

2i<>  (t  r  -  T  ) 

“  2q  -  —7 — (non-melt  surface)  (76) 

(iror  Iq/U ) 

The  reader  is  cautioned  that  non-melting  sliders  undergoing  extremely 
large  sliding  distances  should  be  modeled  as  a  steady-state  conduction 
problem  governed  by  convection  heat  transfer  coefficients  at  the 
boundary. 

Once  it  is  determined  which  surface  melts,  equations  (69)  and  (70,) 
are  equated  to  the  proper  flux,  equations  (74)  and  (76).  The  resulting 
two  equations  can  be  solved  for  film  thickness  h  and  surface  tempera-' 
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ture  Tf  or  Ts,  depending  upon  which  one  has  been  replaced  by  the  melt 
temperature  Tra.  The  'results  show  that  the  same  governing  equations 


apply  to  a  melting  slider  or  a  melting  fixed  plane: 

(Ts,f  ~  T0)  -  [(Tn  -  T0)  ♦  4>3]/(  1  ♦  b2h)  C7TJ 

bjb2h5  +  b^h^  +  b2[(Tm  -  T0)  -  <p3 ] ia  -  2$j  *  0  (73) 

where 

b,  »  o[ (H/c)  +  ( Tm  -  T0)]/ayi26  (79) 

b2  -  2k'(ira'lc/U)_1/2/k  (80) 

$3  *  yU2/2k  (81 ) 


The  solution  to  the  problem  can  be  simplified  if  the  non-melting 
surface  is  conducting  or  non-conducting.  The  total  heat  flux  from  the 
film  is  divided  between  the  melting  and  the  rion-melting  surfaces: 


93  ♦  Qf  -  qm  +  9nm 
From  equations  (69)  and  (70), 
q3  ♦  qf  -  yU2/h 
Now  the  heat  transfer  ratio 

qm  p<jh3[H  ♦  c(Tm-T0}]  /  ira'tf!\  1/2 

t.  ‘  \  u  ) 


(82) 

(83) 


(84) 


represents  two  extreme  approaches:  (i)  qn/qnm  <*  1  “here  most  of  the 
generated  heat  is  conducted  into  the  non-melting  surface,  and  (II) 
qm/qnm  >>  1  where  most  of  the  generated  heat  goes  to  melting  the 
opposing  surface.  Although  many  conditional  factors  affect  thp  ratio 
in  equation  (57),  the  author  has  found  that  surfaces  with  high  thermal 
conductivity  (metals)  tend  to  give  a  solution  that  forms  a  low  ratio 
although  high  sliding  velocities  for  all  material  pairs  will  do  the 
same.  Using  the  words  "conducting  surface"  and  "non-conducting 
surface"  to  characterize  the  approaches,  two  extreme  solutions  can  be 
obtained  for  the  film  thickness.  If  qm/qnm  <<  1 •  then 

uU2/h  -  qnm  (851 

For  most  cases  equation  (77)  reduces  to 


(*a,r  '  To)  -  ('4  -T0) 

Thus',  the  film  thickness  for  a  conducting,  non-melting  surface  is 

»  -  (86, 

2k»(Tm  -  T0} 

Similarly,  if  q0/cinm  >>  ?  *  then 

uU2/h  -  qm  (87) 

Thus,  the  film  thickness  for  a  non-conducting,  non-melting  surface  is 

(MW, 2,  '  T/* 

h  -  - rr  (88)  ' 

po[h  ♦  c  (tb-t0)] 

Once  equation  (86)  or  equation  (88)  is  selected  to  find  the  film 
thickness,  equation  (84)  can  be  used  to  verify  that  the  corresponding 
heat  transfer  ratio  is  satisfied.  ■  If  the  ratio  is  near  unity,  the 
fifth  order  polynomial,  equation  (78),  must  be  solved  for  the  film 
thickness. 

3.  FRICTION  FORCE 

Again  neglecting  the  contribution  of  the  pressure  gradient  to  the 
flow  field,  tne  force  per  unit  area  on  the  slider  is  simply 


The  coefficient  of  melt  friction  is 


The  two  extreme  solutions  for  the  film  thickness  yield 


SK’lT*  -  T Qy 


(conducting  surface) 


-  (u?U2p[H  ♦  c(Tm  -  T„); 

f  *; - iw  . 


(non-conducting 

surface) 


The  Importance  of  material  thermal  properties  is  evident. 
Furthermore,  low  sliding  velocities  suggest  an  increasing  coefficient 
of  friction  with  speed,  equation  (92),  whereas  high  sliding  velocities 


suggest  a  decreasing  coefficient  of  friction  with  speed,  equation  (91). 
The;*e  effects  have  been  noted  in  data  for  the  dry  sliding  of  polymers 

C  353 [ 36 ] -  ' 

4 WEAR 

The  material  wear  is  represented  by  the  amount  of  melt  squeezed 
out  from  the  contact  area,  A,  and  is  given  by  equation  (64).  Expressed 
as  the  volume  removed  per  unit  distance  travel 

v(cm3/cm)  »  —3 -  (93) 

5.  EXAMPLES 

The  melt  lubrication  model  is  now  applied  to  two  examples:  (1)  a 
copper  rotating  band  melting  on  a  steel  gun  barrel;  (2)  an  ice  skater. 
Each  example  represents  conducting  surfaces. 

Rotating  bands  have  no  flow  in  the  circumferential  direction  while 
an  ice  skate  is  very  narrow.  In  each  case  8  is  infinity.  The 
non-melting  steel  has  a  characteristic  contact  length  ic  equal  to  the 
rotating  band  width  when  it  is  the  barrel,  but  it  has  a  characteristic 
contact  length  tc  equal  to  the  skater's  stroke  when  It  is  the  slider. 
Table  9  lists  appropriate  variables.  Melted  copper  properties  were 
obtained' from  Sternlicht  and  Apkarian  [ 45] -  A  summary  of  the  two 
examples  is  given  in  Table  5. 

The  Reynolds'  number  R  is  indeed  small,  and  the  temperature  across 
the  film  differs  little  from  the  melt  temperature.  Thus,  assumptions 
(1)  and  (2)  are  Justified.  Density  changes  sufficient  to  cause  a 
thermal  wedge  action  is  remote.  It  remains  to  verify  that  the  pressure 
gradient  contribution  to  the  energy  equatioh  (66)  is  negligible  when 
compared  to  the  Couette  contribution.  The  parameter  in  equation  (67) 
has  a  magnitude  of  approximately  1 0” 3  for  the  two  examples  above. 

Montgomery  [10]  found  indirectly  that  the  friction  coefficient  of 
gilding  metal  (90*  Cu;  10*  Zn)  rotating  bands  in  155  mo  howitzers 
decays  rapidly  to  f  -  0.020  within  12  cm  of  travel.  Using  his  speed, 
9000  cm/sec,  at  which  this  value  is  reached  and  the  band  width,  2  6  cm, 
the  calculated  value  of  the  friction  coefficient  is  0.024.  Montgomery 
notes  that  "the  coefficient  of  friction  slowly  falls  with  increased 
travel  beyond  this  point  but  this  probably  was  an  artifact  caused  by 


neglecting  the  wear  of  the  rotating  bands."  Although  this  explanation 
is  a  factor,  it  is  significant  that  the  expression  for  the  friction 
coefficient,  based  on  equation  (91),  decays  as  the  square  root  of  the 
velocity.  This  phenomenon  would  be  expected  as  the  projectile 
accelerates  down  the  barrel. 


TABLE  4. 

VARIABLES  FOR 

MELT  LUBRICATION 

EXAMPLES 

Melted 

Melted 

Solid 

Parameter 

Copper 

Ice 

Steel 

U  (cm/sec) 

3.0  x  104 

100 

a  (N/cm2) 

3.4  x  101* 

100 

H  (N-cm/kg) 

1  17  x  106 

1.90  x  107 

Tm  -  T0  CC) 

1055 

2 

a  (cm2/sec) 

0.9 

1.30  x  10-3 

0.09 

c  (N-cm/kg-°C) 

4.3  x  104 

4.3  x  105 

k  (N/sec-°C) 

354 

0.55 

33 

p  (kg/cm3) 

8.9  x  10~3 

1  x  10_3 

u  (N-sec/cm2) 

3-3  x  10-7 

1.7  x  10-7 

i  (cm) 

ic  (cm) 

0.7 

0.3 

0.7  (fixed) 

150  (slider) 

bj  (oc/cm1*) 

2.6  x  1014 

4.1  x  10H 

' 

b2  (cm-1) 

72 

183 

<«3  (°C) 

0.42 

3.1  x  10-3 

TABLE  5.  SUMMARY  OF  EXAMPLES 


Example  ,  h  (cm)  Ts>f  -  Tm  (°C)  f  R 

rotating  band  1.1  x  10-5  0.4  0.026  90 

skater  1.7  x  10'” 5  3  x  10“3  0.011  0.1 


SECTION  V 


FRICTION  THEORY  FOR  METAL  PIN-OK-DISK  DEVICES 

1 .  INTRODUCTION 

Friction  theories  have  been  based  on  material  mechanical 
properties  and  limited  to  low  speeds.  However,  in  most  engineering 
applications  frictional  heating  is  significant.  It  has  been  known  for 
decades  that  dry  rubbing  between  bodies  produces  very  high  temperatures 
at  the  asperities  sufficient  to  melt  the  material  [T4].  Early  high 
speed  sliding  experiments  at  the  Franklin  Institute  DGTJ  have  found 
molten  layers  between  metal  pins  rubbed  on  steel  disks,  the  data 
subsequently  summarized  by  Montgomery  [9].  Recent  workers'  C35][38]C39] 
have  reported  evidence  of  melt  layers  when  polymer  pins  are  rubbed  on 
steel  and  glass  disks  (or  cylinders)  where  the  relative  sliding 
velocities  are  on  the  order  of  one  m/s. 

Stiffler  [12]  derived  a  melt  theory  to  explain  friction  and  wear 
between  sliding  surfaces,  applying  it  to  projectile  rotating  bands. 

The  essential  features  of  the  theory  are  the  following. 

1)  Dry  rubbing  at  sufficient  loads  and  speeds  will  produce 
enough  heat  to  melt  one  surface. 

2)  Once  a  layer  of  melt  forms  a  steady-state  condition  exists  in 
which  viscous  friction  supplies  the  heat  to  melt  one  member 
of  the  sliding  pair. 

3)  The  load  support  is  governed  by  classic  squeeze  film 
principles,  but  the  melt  is  continually  supplied  to  maintain 
a  constant ■ film  thickness. 

4)  The  melt  thickness  must  exceed  the  non-melting  surface 
roughness. 

The  purpose  of  this  paper  is  to  extend  the  rectangular  theory  to 
cylindrical  contacts  and  apply  the  results  to  pin-on-dlsk  data.  The 
theory  development  will  parallel  Section  IV. 

2.  ANALYSIS 

Consider  a  round  pin  sliding  on  a  rotating  disk  in  which  the  disk, 
has  a  higher  melting  point  than  the  pin.  The  following  assumptions  are 
employed  to  describe  the  melt. 


1)  The  fluid  is  laminar  and  incompressible. 

2)  The  pressure  is  constant  across  the  thickness  h,  and  the 
density  and  viscosity  are  constant  throughout  the  film. 

3)  The  pin  is  melting  with  resolidification  after  squeezing  from 
the  contact. 

4)  The  surfaces  remain  parallel. 

5)  Film  thickness  and  heat  penetration  into  the  solid  are  ,°“.il 
compared  to  the  pin  size. 

6)  Quasi-steady  state  conditions  are  reached  in  the  film. 
Assumptions  (5)  and  (6)  are  justified  by  the  theoretical  heat  transfer 
work  of  Landau  [100]  and  summarized  ih  Section  IV. 

Although  the  round  pin  suggests  cylindrical  coordinates  and  the 
sliding  velocity  boundary  condition  suggests  rectangular  coordinates,  a 
difficult  two-dimensional  problem  can  be  avoided  by  using  the 
superposition  principle  [102],  i.e.,  the  non-linear  inertia  terms  are 
neglected  in  lubrication  fluid  mechanics.  Thus,  the  pressure  and 
velocity  fields  for  the  "Couette  effect"  due  to  pin  translation  can  be 
added  to  the  fields  for  the  radial  "squeeze  effect." 

The  Couette  effect  develops  no  pressure  field  or  load  support. 

The  velocity  i3  3imply 


where  x  denotes  the  direction  of  disk  sqrface  velocity  U  and  y  denotes 
the  direction  perpendicular  to  the  surfaces. 

The  momentum  equation  for  the  squeeze  effect  is 


Integrating  across  the  melt  film  for  a  stationary  di?k  and  pin, 

'  vr  -  1—  (y2  -  yh)  ,  •  (96) 

2\i  dr 

.The  continuity  equation  is 

1  a(>r  V  .  iiadt  .  o  (97) 

r  3r  3y 

Substituting  equation  (96)  and  integrating  across  the  film  gives  the 
reduced  Reynolds.'  equation  for  parallel  surfaces: 


68 


1  £_ 

r  dr 


R] 


12y  f  ,h 
Ph3  lpv^0 


(98) 


p(R)  »  pa 

dp(0)  „ 

dr 

The  term  Vy  at  the  surface  is  always  interpreted  as  the  rate  of  surface 
approach  3h/3t,  i.e.,  a  squeeze  film.  However,  in  this  case,  the  term 
p Vy  »  til  must  be  interpreted  as  the  rate  of  mass  addition  per  unit  area 
from  the  melting  pin. 

The  solution  is  given  by 


3  um  p  p 

p  "  pa  -  7^3  (r2  '  p2) 

Integrating  over  the  pin  area  for  the  load  capacity, 


(99) 


(100) 


or 


3umR* 

2ph"3 


(101) 


where  a  is  the  average  pressure  between  the  pin  and  disk. 

It  is  viscous  dissipation  in  the  melt  layer  that  supplies  the 
energy  to  continually  generate  pin  melt.  In  addition,  energy  is  lost 
from  the  film  by  way  of  transient  heat  conduction  into  the  non-melting 
disk  surface.  The  film  thickness  (degree  of  melt)  depends  on  the 
distribution  of  heat  between  the  pin  and  disk.  Since  the 
characteristic  time  of  the  film  heating  (h2/a)  is  much  smaller  than  the 
pin  contact  time,  the  film  heating  is  quasi-static.  Furthermore,  only 
gradients  perpendicular  to  the  surfaces  are  important  (assumption  5). 
Thus,  the  energy  equation  for  the  melt  layer  reduces  to 


d2T 


dy 


-[e)Mr)2]-° 


(102) 


•r.V.  f.  A  •*:  r;;*.  v.  v/.v 
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Now  in  all  practical  cases  the  pin  translational  velocity  vx  is  much 
greater  than  the  squeeze  velocity  vr,  and  viscous  heating  from  the 
latter  is  neglected.  Assuming 

dVp/dy  dh2<j 
— C — L  =  -  «  ! 

dvx/dy  iruRU 

the  solution  to  the  energy  equation  is 

ulj2 

T  -  Tra  -  (Td  '  Tm)(y/h)  -  (y2  -  yh)  (103) 

where  Td  and  Tm  are  the  surface  temperatures  of  the  disk  and  melting 
pin  respectively. 

The  heat  transfer  per  unit  area  out  of  the  film  is  given  by  the 
temperature  gradient  at  the  surface: 

qp(pin)  -  +k[dT/dy ]0 

qd(disk)  »  -kfdT/dy]^ 

or 

qp  -  -k(Tm  '  Td]/h  +  pU2/2h,  (104) 

Qd  -  +k(Tm  "  Td)/h  +  pU2/2h  (105)' 

What  is  the  heat  transfer  in  equations  (104)  and  (105)?  For  the 
melting  pin  the  work  of  Landau  [100]  is  important.  He  found  that  the 
problem  of  a  heat  flux  imposed  on  a  surface  where  the  melt  is 
continually  removed  will  reach  a  steady  state  solution  within  a  time 
scale  much  less  than  the  contact 1  time.  The  state  is  characterized  by  a 
steady  inward  migration  of  a  fixed  temperature  distribution  within  the 
solid,  and  the  heat  transfer  is  given  by 


dp  -  m[H  ♦  c(Tm  -  T0j]  '  (106) 

At  the  disk  surface  the  heat  transfer  for  transient  conduction  into  a 
semi-infinite  solid  is 


k’(Td  “  T0) 


(107) 


\ 
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where  the  prime  represents  disk  properties.  For  a  pin  with 
characteristic  length  £0  sliding  on  a  disk  at  relative  velocity  U,  the 
characteristic  contact  time  tc  -  lc/ U.  Using  the  average  heat  flux  for 
the  contact  time. 


<Jd 


2k’ (Td  -  T0) 
(iraUo/U)1 


(108) 


Since  contact  length  varies  across  the  pin  an  average  contact  length  is 
defined: 


2R 

ir/2 


■  it/ 2 

CO30d9 

0 


4 

— ,  R 

7T 


(109) 


Equations  (104)  and  (105)  may  be  solved  now  for  the  two  unknowns:  disk 
surface  temperature  Td  and  film  thickness  h; 


(Td  -  T0)  -  [(Tm  -  T0)  ♦  *3]/(l  *  <f»2h) 

(110) 

<h4>2h5  +  h4  +  <f»2((Tm  “  T0)  -  ij>3)h  -  2*3  «  0 

(111) 

where 

*1  -  2<j[(h/c,)  +  (Tm  -  T0)]/3ohjR2 

(112) 

<t>2  -  2k’(ifa’lc/u)"1/2/k 

(113) 

*3  -  yU2/2k  ■ 

(114) 

For  all  practical  cases, 

(Td  “  T0)  -  (Tm  -  T0) 

(115) 

Equation  (111)  for  the  film  thickness  has  a  simple  solution  for 

two  extreme  conditions:  (1)  qp/qd  <<  1  where  most  of  the  generated  . 
heat  is  conducted  into  the  disk,  called  a  conducting  disk  surface;  (2) 
Qp/Qd  >>  1  where  most  of  the  generated  heat  goes  to  pin  melt,  called  a 
non-conducting  disk  surface.  This  ratio  is 

q„  p.h3(H  .  c(T„  -  To)] 

H  3»nVfTn  "  T»)  '  0  ' 
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While  there  are  many  parameters,  the  film  thickness  i3  the  predominant 
term.  These  ^iim  thicknesses  are  very  small  and  the  ratio  is  much  less 
than  one  for  all  pairs  examined  by  the  author.  For  example,  low 
melting  nylon  pins  sliding  on  a  glass  disk  at  one  m/s  gives  a  ratio  of 
1CT2. 

From  equations  (104)  and  (105) 

qp  ♦  qd  =  uU2/h  (117) 

For  a  conducting  disk,  qp/q^  <<  1,  the  film  thickness  is 


uU2(q'R/U)1/2 

k'(Tm  -  T0 I 


(118) 


For  a  non-conducting  disk,  qp/qq  >>  1,  the  film  thickness  is 

1/4 


h 


3 ( uUR ) ‘ 


2po[H  *  c(Tm  -  T0 J  J 

3.  FRICTION  FORCE 

The  friction  force  per  unit  area  is 

pU 


(119) 


(120) 


The  coefficient  of  melt  friction  is 


f  -  & 
ho 


The  two  extreme  conditions  for  the  film  thickness  yield 
,  *'(*«  -  To) 


o(a’RU)1/2 


2p2U2p (h  +  c(Tm  -  T0) 
3R2a3 


(conducting  disk) 
1/4 


(non¬ 
conducting  disk) 


(121) 


(122)  ■ 


(123) 


4.  WEAR 

The  material  wear  is  represented' by  the  amount  of  melt  ill  squeezed 
from  the  contact  area  and  is  given  by  equation  (101).  Expressed  as  the 
volume  removed  per  unit  distance  of’ travel: 

2noh3 


V(cm3/cm) 


3pU 


(124) 
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5.  RESULTS 

Perhaps  the  most  comprehens i ve  experimental  study  of  metal 
friction  and  wear  at  high  sliding  speeds  was  done  at  the  Franklin 
Institute  from  1946  to  1956.  A  pin-on-disk  device  was  used  to  slide 
metal  pins  aganst  a  spinning  steel  disk  at  speeds  up'  to  548  m/s.  The 
pins  were  primarily  copper  and  gilding  metal  (90  CU;  10  ZN)  although 
several  runs  exist  for  aluminum,  zinc,  and  nylon.  A  radial  pin 
displacement  on  the  disk  assured  that  all  data  were  single  pass.  This 
fact  is  important  for  comparisons  between  theory  and  experiment. 

A  calculation  of  the  heat  transfer  ratio  showed  that  the  disk  is 
conducting  for  all  test  parameters.  Therefore,  a  normalized  coefficient 
of  friction,  based  on  equation  (29),  is  defined; 


f 


f0q(q'RU)1/2 
k,(Tm.-  T0) 


Cl  25) 


where  fQ  is  the  measured  coefficient  of  friction.  The  data  plotted 
against  f  should  reduce  to  one.  Thermal  properties  [103]  of  the  gun 
steel  disk  can  vary  significantly  with  temperature: 


T(°C) 

27 

250 

400 

600 

1250 

,  N  ,  „  „  B'tu 

k'( - )  -  1.73  x - 

sec°C  hr  ft  °F 

43 

40 

36 

33 

32 

a’(cm2/sec) 

0.122 

0.092 

0.073 

0.058 

0.063 

Since  the  heat  flux  is  defined  at  the  disk  surface,  the  thermal 
conductivity  in  equation  (125)  should  be  the  value  at  the  pin  melting 
temperature'.  However,  the  disk  diffusivity  covers  the  range  from  Tm  to 
T0 ,  and  an  average  diffusivity  (not  the  diffusivity  at  the  average 
temperature)  is  used.  It  is  assumed  that  the  disk  bulk  temperature  T0 

is  27 °C.  - 

A  plot  of  normalized  coefficient  of  friction  versus  pin  pressure 
is  shown  in  Figure  27  for  a  wide  range  of  sliding  velocities.  The 
majority  of  the  data  represents  copper  pins  with  a  diameter  of  u.2031 . 
cm-  The  available  data  for  aluminum,  zinc,  and  nylon  was  limited  to 
274  m/sec  with  a  pin  diameter  of  0.3556  cm.  There  is  a  rather 
remarkable  fit  for  copper  at  the  higher  pin  pressures  and  a  tight 
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Figure  27.  Experimental  over  Theoretical  Coefficient  of  Friction  as  a  Function  of 
Sliding  Speed  and  Load  per  Unit  Area 


correlation  with  sliding  velocity  at  all  pin  pressures.  The  measured 
coefficient  of  friction  at  lower  pin  pressures  is  smaller  than  the 
coefficient  predicted  by-  theory.  The  reason  is  the  disk  roughness. 

Shugarts  and  Rippei  [101]  reported  that  the  Franklin  Institute 
test  disk  had  an  rms  roughness  of  5  -  13ucra  (2  -  in  the  direction 

of  travel.  The  copper  film  thickness  was  calculated  from  equation 
(118),  using  a  melt  viscosity  of  3  3  x  10"?  N-see/co^;  h  -  8|icm.  Under 
very  high  pin  pressures  the  disk  asperities  are  greatly  deformed  and 
the  surface  roughness  is  reduced.  At  lower  pin  pressures  the  copper 
film  thickness  and  the  disk  asperities  are  comparable  and  the  film  is 
no  longer  uniform  ever  the  pin  surface.  Thus,  the  pressure  is  larger 
than  t.ne  calculated  value  using  the  pin  area,  and  the  coefficient  of 
friction  will  be  smaller.  In  effect  the  normalized  coefficient  of 
friction  for  copper  at  low  pressures  is  less  than  the  theoretical  value 
because  the  melt  layer  i3  incomplete,  and  the  pressure  a  is  based  on 
the  pin  area,  not  on  the  actual  contact  area.  The  lower  melting  point 
nylon  and  aluminum  produce  thicker  melt  films  and  give  satisfactory 
results  at  the  lower  pin  pressures.  However,  the  data  for  zinc  is 
inconsistent  witn  these  findings. 

The  wear  rates  of  the  copper  pins  were  tabulated  also  for  the 
Franklin  Institute  tests  [93.  These  experimental  wear  rates  were 
approximately  10"®  cm3/cm.  When  the  theoretical  wear  rate  from 
equation  (31)  was  calculated,  it  was  found  to  be  approximately  10"® 
cm3/cm.  Wear  rates  are  difficult'  to  compare  because  of  the  sensitivity 
of  the  wear  rate  to  the  melt  viscosity.  It  is  certain  that  the 
viscosity  is  much  higher  than  the  number  used  here  because  of  the  very 
high  pin  pressures.  Also,  the  high  shear  rates  are  a  factor.  Mo 
information  appears  to  be  available  on  how  these  two,  parameters  affect 
the  viscosity  of  metals. 

Contact  Area ■ 

The  friction  coefficient  can  be  lowered  by  increasing  the  contact 
stress  either  by  (1)  decreasing  the  pin  area,  or  (2)  increasing  the 
number  of  contacts,  i.e.,  a  surface  with  a  roughness  exceeding  the  melt 
thickness.  Assume  N  contacts  of  radius  a.  Then  equation  (122) 
becomes 


(126) 


Nira3/2k’(Tm  -  T0) 
w(a'U)1/i; 


where  W  is  the  known  load.  The  coefficient  of  friction  decreases  with 
the  number  and  size  of  the  contacts  until  the  flow  pressure  <jf  of  the 
softer  surface  is  reached.  The  latter  state  is  achieved  for  most  rough 

surfaces  under  light  loads  and  slow  sliding  spe-'d.n  where  the  melt 

thickness  is  much  less  than  the  roughness  height.  Gupta  and  Cook  [104] 

have  investigated  several • pairs  of  mating  rough  surfaces  and  found  that 

the  number  of  contacts  is  nearly  proportional  to  the  load  while  the 
size  of  the  contact  is  nearly  independent  of  the  load,  assuming  round 
asperiti  of  equal  size.  They  found  that  the  ratio  of  : low  pressure 
to  hardness , falls  between  0.50  -  0.55.  Replacing. the  contact  stress 
with  the  flow  pressure  in  equation  (122), 


where 


f 


k’(Tm  -  T0) 
Of (a'aU ) ' 


(127) 


a  .  4R  x  10“3//c  (128) 

The  topographic  index  c  introduced  by  Gupta  and  Cook  characterizes  the 
roughness:. 


E’ 

C  -  - 

8Hv/2Jf 


(129) 


where  0  is  the  asperity  radius  of  curvature  and  N  is  the  number  of 
peaks  prior  to  .loading.  Thus,  the  coefficient  of  friction  is  a  minimum 
when  the  contact  stress  approaches  the  pin  flow  pressure,  and  the  size 
of  the  individual  contact  spots  is  as  large  as  possible  while 
maintaining  the  high  contact  pressure.  In  theory  frictional  behavior 
can  be  improved  by  dimpling  the  non-melting  surface. 

Tanaka  and  Uchiyama  [39]  udied  the  effects  of  sliding  speed  on 
the  friction  and  wear  of  c  /stalline  polymers.  A  pin-on-disk  device 
was  used  for  speeds  over  the  range:  10  -  300  cm/s.  Overwhelming 
evidence  was  presented  to  support  a  melt  phenomena.  Melt  depth  of 
polymer,  pins  on  a  glass  disk  was  on  the  order  of  1.0“2  cm  while  the 
depth  with  a  steel  disk  was  only  several  microns.  The  heat  ratio  in 
equation  (23)  was  much  less  than  one  for  both  glass  and  steel  disks.  ' 


Thus,  the  theoretical  coefficient  of  friction  is  given  by  equation 
(122).  Table  6  lists  the  polymers  and  compares  the  theoretical  and 
experimental  coefficients  of  friction  for  a  glass  disk.  Glass 
properties  [105]  are  k'  =■  1.17  N/sec-°C  and  a’  =  0.006  cra^/sec. 

It  is  observed  that  the  theoretical  coefficient  of  friction 
exceeds  the  experimental  values  by  a  factor  of  three.  The  measured 
coefficients  are  lower  than  predicted  because  all  tests  were  multiple 
pass.  Tanaka  and  Uchiyama  reported  polymer  melt  transfer  to  the  disk 
where  polyacetal  was  the  worse  offender.  They  obtained  a  photograph  of 
the  contact  surface  through  the  glass  disk.  .  There  existed  a  number  of 
narrow  striations  parallel  to  the  rubbing  direction  on  the  frictional 
surface.  These  bands  were  identified  as  non-contact  areas 


TABLE  6.  FRICTION  COEFFICIENTS  AND  WEAR  RATES  FOR 
SEVERAL  POLYMERS  AGAINST  GLASS.  (Experimental  data 
from  Tanaka  and  Uchiyama  [39].) 


Polymer 

M.C.  ' 
(°C) 

f 

(exp, 

f 

.)  (th.) 

f 

V  (exp.) 

V  <th. ). 
cm/ cm 

Low  Density 
Polyethylene* 

115 

0.9 

2. 

3 

0  39 

4  x  10"6 

— 

High  Density  . 
Polyethylene** 

135 

0.7 

2. 

1 

0.33 

4  x  10-8 

5  x  10'5 

Polypropylene* 

167 

1.3  ' 

3. 

I 

5 

0.37 

i  x  ra~7 

3  x  10~5 

Nylon**  . 

250 

1.3 

4 

3  ' 

,0  30 

1  x  10"6 

1  ‘X  IQ-*1 

Polyacetal** 

1 75 

0.6 

2, 

9 

0.21 

3  x  10-7 

2  x  10"5 

*120  cm/s 
»*  250  cm/s 

due  to  polymer 

lumps  and 

film 

which 

had 

adhered 

to  the  disk 

Therefore,  full  contact  does  not  occur,  and  the  actual  contact  pressure 


culated  from  the  pin  size.  The 
t  area  as  less  than  one  half  of  the 


is  much  higher  than  the  pressure  cal 
author  visually  estimates  the  contac 
pin  cross-sectional  area.  The  unknown  disk  bulk  temperature  could  be 
Important,  especially  for  multiple  pass,  tests.  . 

Theoretical  coefficients  of, friction  for  the  steel  disk  are 
approximately  15  -  30.  These  unusually  high  values  are  a  sign  that  (1) 
very  thin  films  are  disrupted  by  the  surface  roughness,  and  high  local 


pressures  exist  in  the  film  or  (2)  low  speeds  and  Ibads  do  not  supply 
sufficient  energy  to  bring  the  contact  into  the  thermal  melt  regime. 
Consider  again  equation  (122).  If  both  sides  are  multiplied  by  the 
sliding  velocity  and  rearranged, 

foU  =  qd  '  (130) 

i.e.,  the  frictional  work  input  per  unit  contact  area  is  equal  to  the 
transient  heat  flux  to  the  disk  surface.  The  contact  stress  must  be 
determined  accurately,  and  the  contact  must  reach  melt  temperature 
before  equation  (112)  is  valid.  .  . 

Table  6  compares  the  theoretical  wear  rate  for  several  polymers 
with  the  experimental  values  from  Tanaka  and  Uchiyama.  Equation’  (124) 
has  been  adjusted  for  the  higher  contact  pressures  and  expressed  as  a 
change  of  pin  length  per  unit  length  of  travel.  The  measured  wear 
rates  are  smaller  than  theory  by  a  factor  of  100.  Polymer 
viscosities  were  determined  by  measuring  the  solidified  melt  depth  on 
the  pins  and  by  using  an  expression  identical  to  equation  (121).  The 
calculated  viscosity  agreed  with  viscometry  within  shear  rate  limits. 
One  possible  explanation  for  the  difference  is  the  measurement  method. 

A  linear  differential  transformer  was  used  to  measure  the  change  In  the 
pin  length.  There  w  .s  no  discussion  of  any  adjustment  made  for  the 
transferred  polymer  lumps  on.  the  glass  disk. 

6.  SUMMARY 

The  proposed  melt  theory  represents  a  new  approach  to  friction  in 
the  thermal  regime.  A  tentative  evaluation  of  the  theory  can  be  made 
based  on  the  limited  studies.  The  theory  predicts  the  coefficient  of 
friction  for  metal  pairs  sliding  at  high  loads  and  speeds  under  single 
pass  conditions.  At  typical  engineering  load3,  surfaces  with  high 
thermal  conductivity,  such  as  metals,  produce  very  thin  melts,  on  the 
order  of  the  surface  rcughne?,3.  Then  the  size  and  pressure  of  the 
contact  i3  uncertain  and  the  measured  coefficients  of  friction  always 
are  less  than  the  calculated  values.  A  comparison  of  experimental  wear 
rates  with  a  melt  theory  of  wear  was  Inconclusive.  Data  on  the 
viscosity  of  metal  melts  at  high  contact  pressures  and/or  rates  of 


shear  are  generally  not  available.  In  the  case  of  multiple  pass 
contact  the  effect  of  transfer  debris  on  idealized  flow  passages  is 
difficult  to  assess. 

In  the  case  of  polymer  pin  friction  and  wear  a  comparison  of 
theory  with  one  particular  experiment  [39]  was  inconclusive.  Multiple 
passes  lead  to  the  presence  of  transfer  (wear)  debris  which  altered  the 
apparent  contact  pressure,  a  necessary  parameter  In  the  theory. 

However,  this  problem  should  not  occur  for  the  single  pass  polymer 
rotating  bands.  More  important,  the  non-linear  behavior  of  polymer 
melts  calls  into  question  an  application  of  Newtonian  flow  for  these 
materials.  The  subject  of  non-Newtonian  friction  theory  will  now  be 
addressed. 
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SECTION  VI 


FRICTION  THEORY  FOR  PLASTIC  ROTATING  BANDS 

In  the  previous  two  sections  melt  friction  theory  was  derived  for 
Newtonian  fluids  (in  particular  metal  melts)  and  applied  to  copper 
rotating  bands  and  metal  pin3.  Good  agreement  between  theory  and 
experiment  was  achieved  to  confirm  the  basic  approach  to  the  problem. 

In  this  section  melt  friction  theory  will  be  derived  for  polymer  melts 
which  display  the  nor.-Newtonian  fluid  behavior  discussed  in  Section 
III.  The  fluid  is  characterized  by  viscosity,  elastic  shear  modulus, 
and  limiting  shear  stress,  all  of  which  can  vary  with  pressure  arid 
temperature.  - 

Although  plastic  rotating  bands  are  circular  the  melt  films  are  so 
thin  that  a  Cartesian  coordinate  system  can  be  used  to  analyze- the 
problem.  Imagine  the  thin  film  cut  and  spread  out  in  the  z-direction 
with  the  x-direction  of  motion  as  the  length.  Since  the  fluid  can  not 
flow  in  the  z-direction  (the  circular  band  is  continuous),  the  analysis 
is  one-dimensional  in  the  plane  of  the  contact.  The  rotating  band  is 
referred  to  as  the  "moving  surface"  while  the  barrel  is  referred  to  as 
the  fixed  surface. 

1 .  MOMENTUM 

The  analysis  is  begun  with  the  usual  five  assumptions  at  the 
beginning  of  Section  IV  which  are  characteristic  of  thin  film  fluid 
mechanics  in  lubrication  theory.  Since  It  is  no  longer  possible  to 
start  with  derived  forms  of  Reyndlds'  equation  for  pressure 
distribution  such  as  equation  (60),  it  is  necessary  to  follow  those 
derivations  for  the  non-Newtonian  lubricant.' 

Neglecting  fluid  inertia  and  body  forces,  the  general  momentum 
equations  are  [46] 

i£  .  iiXX  *  llXi  (131). 

3x  3y  3z 

iE  .  ilS* llXJL  (132) 

3y  3x  3z 
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033) 


_3p  m  _9jjrz 

3z  3x  3y 

Since  the  film  is  thin  in  the  y-direction,  3p/3y  =  0,  Al3d, the  film 
plane  is  one-dimensional  in  the  x-direction  or  3p/3z  =  0.  In  addition, 
there  are  no  shear  stresses  in  the  z-direction.  Thus,  the  above  ■ 


reduces  to  (xXy  »  i) 

3p  3t 
3x  '  3y 


(13*) 


a 

0  *  — 
3x 


0  35) 


With  the  exception  of  Conry  [106]  equation  035)  is  usually  neglected; 
the  consequences  are  major  as  discussed  below. 

Integrating  with  respect  to  y  and  assigning  xj  a3  the  shear  stress 
at  the  fixed  surface', 


dp 

^  T1 


0  36) 


Now  the  Urge  contact  pressures  and  large  relative  sliding  at  the 
contacts  will  assure  that  the  limiting  shear  stress  x^  will  occur  at 
the  sliding  surface  boundary: 


T1  * 


037) 


The  lubricant  shear  behavior  is  described  above  by  the  non-linear 


Maxwell  model: 


du  U  dx 
dy  G  dx 


*  5 


038)- 


From  equation  035)  it  would  appear  that  the  elastic  term. is 
inoperative.  However,  let  us  examine  the  elastic  term  more  closely, 
.If  a  small  shear  strain  c'is  considered  to  bq  made  up  of, a  viscous 
component  cv  and  an  elastic  component  ce>  we  have 

C  -  Cv  * 

Then  from  Newton's  law  of  viscosity  and  Hooke's  law  of  elasticity,  the 
rate  of'  strain  is 
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pressure  distribution.  Employing  the  integral  form  of  the  continuity 
for  one-dimensional  flow  with  m  as  the  mass  rate  of  melt  per  unit 
area: 


(147) 
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Jo 

and  equation 

(146),  • 
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/I  Uk  dp  \ 
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dx  [ 

5  u  G  dx  / 

'  3 

dx  u  2  /J 

(143) 


Define  dimensionless  variables: 
p  -  pBL/W 
x  »  x/L 

where  S  i3  the  rotating  band  circumference ,  L  is  the  rotating  band 
length,  and  W  is  the  load.  After  employing  average  values  (defined  at 
the  average  stress)  for  viscosity  y,  shear  modulus  G,  and  limiting 
shear  stress 

•T-  f  ( 1  -  Bi  Q  )(  B2  -  1  )  I  *  -  B0  (149) 
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/  „  dP  ,  \  I 

1 

1  where 

dx 

t 

ca 

Nil 

{B2  t*  -  vj 

B0  2Lmy/pTL1h‘- 
B,  *  IkyW/GBL2 

B2  =>  2hW/3iLBL2  . 

Boundary  conditions  for  equation  (149)  are  the  foliowlnjgf 
p(0)  =*  0  ' 

p(1)  -  0 

A  propellant  gas  pressure  exists  on  the  front  side  of  the  band,  but 
the  pressure  i3  symmetrical  and  will  not  contribute  to  the  load 
support. 

Integrate  the  above  equation . once  to  arrive  at 


1  -  B, 


)  (b2 


-  i  -  -  Bnx  ♦  Ci 


(150) 


The  pressure  gradient  can  be  found  by  using  the  quadratic  formula.  A 
second  integration  gives  the  pressure  distribution.  However,  the  first 
integration  constant  is  embedded  within  a  radical  which  makes  it 
impossible  to  apply  boundary  conditions  and  obtain  a  complete 
analytical  solution.  Numnerical  methods  must  be  employed'.  On  the 
otherhand  a  simplification  can  be  carried  out  which  should  be  valid  for 
all  the  physical  cases  which  will  occur  in  practice.  Rearranging 
equation  (150): 


(«gf  +  /!_  «.  L_\  ^  .  . 

Vdx/  \B;  B2  '  dx 


x  +  C2 


(151) 


The  pressure  gradient  has  been  normalized  so  that  its  magnitude  is  on 
the  order  of  one;  thus,  the  constant  parameters  dictate  the  order  of 
the  two  terms  on  the  left.  If  either  Bi  «  1  or  B2  «  1,  both  of  which 
are  probable  for  melt  friction  cases,  the  first  term  can  be  neglected 


dp  Bn 

dx  “  (Bi+B2)  X  *  C3 


(152) 


Integrating 


Applying  boundary  conditions, 


—  +  C3x  ♦  C4 


(153) 


2( B; +B2) 


[X2  -  x] 


Now  the  rotating,  band  must  carry  load  W.  Therefore, 


f1  -  - 

p  dx  -  1 
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Completing  equation  (155), 
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The  above  equation  relates  the  load  to  the  mass  rate  of  melt  (and  film 
thickness).  The  energy  equation  to  follow  will  give  a  relationship 
between  the  mass  rate  of  melt  and  film  thickness.  Together,  the  film 
thickness  can  be  found. 

2.  ENERGY 

The  energy  equation  for  an  element  of  fluid  is  a  balance  between 
heat  conducted  into  the  element,  the  heat  generated  through  viscous 
forces,  and  the  rate  of  change  of  stored  energy: 


DT  . 

pc  —  =  k 
Dt 


32T  32T 

3x2  *  3y2 


+  tF(t) 


(157) 


where  the  latter  term  is  the  viscous  dissipation  for  non-Newtonian 
fluids  [87].  By  dimensional  reasoning  the  sliding  process  is  quasi 
steady-state  if  the  characteristic  time  of  contact  is  much  greater  than 
h^,+  diffusivity.  Using  the  assumption  [3^]  that  conduction  is  much 
greater  across  the  film  than  in  the  direction  of  motion  x, 

32T 


pcu 


^T 

3x 


k  +  tF(t) 
3y<- 


(158) 


For  the  non-linear  Maxwell  model  of  the  fluid  the  term  F(t) 
represents  the  viscous  strain  rate.  Since  a  model  has  been  adopted 
that  displays  a  Newtonian  viscous  strain  rate  up  to  the  point  of  slip 
at  the  moving  surfce  wnere  the  shear  stress  is  limited  to  t^, 

F(t)  -  t/p  (159) 

In  addition,  the  shear  stress  and  fluid  velocity  expressions  have  been 
derived  previously: 
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Introducing  dimensionless  numbers 
x  -  x/L 
y  »  y/h 
p  -  pBL/W 
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X  *=  x/xL  . 

the  above  equations  are  written  as 

p.„  s  5 .  |!l .  72  ii-2”2 
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du  hx 
dy  “  uU 


:it  (l  -  Bj  — P-  )  x 
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(163) 
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where 


.  .  h  Uhpc 

P’p  (reduced  Peclet  No.)  =  —  - 

b  k 

The  general  order  of  magnitude  of  tne  terms  in'  indicate  that  it 
will  be  much  less  than  one  while  the  very  thin  melt  thicknesses 
compared  to  the  rotating  band  size  also  indicate  that  the  reduced 
Peclet  number  and  B2  are  much  less  than  one.  Then 

x  s  1  (165) 


du 

dy 


hx 
=  uU 

d2T  xj2h2 

— S*.  +  — « - 

dy^  yk 


(166) 


(167) 


A  small  Peclet  number  not  only  reduces  the  partial  differential 
equation  to  an  ordinary  differential  equation  but  also  uncouples  the 

velocity  functional  from. the  energy  equation.  _ _ 

Let  Tf  be  the  film  temperature  at  fixed  surface1  (barrel)  and  Ta  be 
the  film  temperature  at  the  sliding  surface  (band).  Film  temperatures 
are  not  expected  to  vary  much  across  the  film;  thus,  parameters  are 
constant  with  respect  to  y.  Integrating  equation  (167)  and  applying 
boundary  conditions  (in  the  form  of  surface  temperatures). 


T  -  Tf  -  (,T3  -  Tf  )y 


068) 


Heat  transfer  per  unit  area  out  of  the  film  is  expressed  in  terms 
of  the  film  temperature  gradient  at  the  surface: 
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qs(x)  =  -k(dT/dy)^  (169) 

qr(x)  *  +k(dT/dy)0  (170) 

or 

Qs  -  -  £  UTs  -  Tf)  -  TL2h2/2pk]  (171) 

Qf  =  +  £  [(Tg  -  Tf)  ■+  TL2h2/2yk]  (172) 


The  boundary  conditions  are  completed  once  the  heat  flux  to  the 
surfaces  is  known.  Heat  is  transferred  by  conduction  into  the  solids 
in  amounts  q^  at  the  barrel  and  at  the  melting  band.  When  slip 
occurs'  at  the  moving  surface,  the  energy  input  i3  tlU  while  the  film 
receives  only  t^Ui  where  is  the  surface  velocity  of  the  fluid.  The 
difference  represents  a  heat  source  due  to  slip.  Thus, 

Qs  +  tl('J  -  Ut  )  -  qM  ’  (173) 

and 

9f  *  9nm  (17*0 

Slip  veioc'ties  can  be  found  by  integrating  equation  (167): 

1  du  _ 

■  Ui  =  U  -=  dy  (175) 

Jo  dy 

or 
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From  Section  IV, 


qm  *»  m[H  +  ,c(Tm  -  T0)j 


(1 77) 
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2k' (Tf  -  T0) 
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(178) 


Since  fs  5  Tm  (melting  slider)  equations  (173). and  (174)  can  be  solved 
for  the  film  thickness  and, barrel  surface  temperature  Tf.  Making  the 
proper  substitutions  into  equations  (1.73)  and  (174)  and  using  equation 
(156)  to  define  the  melt  rate  m, 


-(Tm  -  Tf)  -  A0h2  *  Ai<h  =  ATh3  +  A2h4  (179) 
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Solving  for  Tf  and  h. 

•** 

i 

{ 

(Tf  -  T0)  -  [(TB  -  T0)  *  Aoh2]/(l •  ♦  A3h) 

(181) 

' 

A2A3h^  ♦  ( A-j A3  +  A2)h3  ♦  (Ai  ♦  AoA3)b2 
-  A3Ai,h  +  A3(Tm  -  T0)  -  A4  -  0 


(182) 


The  solution  to  the  problem  can  be  simplified  if  the  surfaces  are 
conducting  or  non-conducting  as  defined  in  Paction  IV.  The  ratio  of 
heat  transfer  to  the  surfaces  is 


SfflL.  ,  h2Ul  ♦  *2h) 
°.nm  A3(  Tf  T0 ) 


(133) 


If  the  ratio  qm/qnn)  «  1t  most  of  the  generated  neat  given  by  r^’J  goes 
to  conduction  into  the  non-melting  surface  (said  to  be  a  conducting 
barrel  surface).  Then 

<Jna  -  Utl  (184) 

Setting  Ai  -  Ao  m  0, 

A0A3h^  -  AjAijh  ♦  A3(Tra  -  T0)  -  Aj,  -  0  (185) 

and  the  quadratic  can  be  solved  directly. 

If  the  rat-o  qnj/qnm  »  1.  most  of  the  generated  heat  goes  to 
melting  the  rotating  band  (said  to  be  a  non-conducting  barrel  surface). 
Then  _ _ J 

dm  ■  Utl  (186) , 
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Substituting  from  equation  (177), 


A2h3  *  Ath2  -  UrL/k  (187) 

3.  FRICTION  FORCE 

Tne  friction  force  on  the  rotating  band  is  the  surface  shear 
stress  multiplied  by  the  contact  area.  For  the  non-Newtonian  polymers 
the  shear  stress  is  simply  the  limiting  shear  stress  T£..  Then  the 
coefficient  of  friction  is 


The  coefficient  of  friction  decreases  fcr  increasing  contact  pressure 
as  with  the  metal  rotating  bands;  however,  the  limiting  shear  stress  is 
both  pressure  and  temperature  sensitive.  No  distinction  is  made 
between  the  "conducting"  and  "non-conducting"  barrels. 

«.  WEAR 

The  material  wear  is  given  by  the  amount  of  melt  squeezed  from  the 
contact  area,  m,  and  is  given  by  equation  (156).  Expressed  as  the 
volume  removed  per  unit  distance  travel: 


( 1 89) 
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SECTION  VII 

EXAMPLE  FOR  POLYMER  BAND 

t.  ROTATING  BAND  (30  PROJECTILE) 

Material:  Polyethcrsul phone 
Length:  L  -  0.5  in 
Circumference:  B  -  37  in. 

Speed:  10**  in/sec 

2.  BARREL  PROPERTIES 

k*  -  4.6  x  tO"4  Btu/sec-in-°F 
a'  *  0.023  in2/see 
Length:  60  in 

3-  BAND  PROPERTIES 

p  -  0.035  lbm/in3 
c  »  0.5  Btu/lbm-°F 
H  -  54  Btu/lbm 
Tm  -  T0  -  470°F 
k  -  9.03  x  10“6  Btu/sec-in-°F 
G  -  15  lbf/ln2  (?) 

•c  -  1  x  10~5  in2/lbf  (?)  ! 

.9  x  10“**  in2/lbf) 
reference  39)  «  10"2 


lyners  is  difficult  to 
licate  the  task.  Based 
e  coefficient  of 


B  (equation  50)  -  16,000  1/°K 
a  (equation  52)  -  0.276  x  10“?  n2/N  ( 
Wo  (estimated  from  similar  polymers  1 
lbf-sec/in2 

w  (equation, 51 )  -  0.45  lbf-sec/in2 
4.  FRICTION 

The  coefficient  of  friction  is  defined  as 
f  -  tl3L/W 

However,  data  on  the  limiting  shear  stress  of  p<i 
find  (if  it  exists).  High  pressure  effects  coaqi 
upon  the  author's  work  [8],  it  was  found  that  tit 
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friction  for  polyethersulphone  rotating  bands  under  conditions  similar 
to  actual  firings  was  f  -  0.016  at  a  stress  of  W/BL  »  20,300  lbf/in2. 
Thus, 

tl  -  325  lbf/in2 

5 .  CONSTANTS 

Aq  *  1.4  x  106  °F/in2 
A}  -  5.6  x  1011  °F/in2 
A2  -  1.0  x  1011  °F/in4 
A3  -  5.3  x  104  1/in 
At)  *  3.8  x  10?  °F/in 

6.  MELT  THICKNESS 

h  (equation  1 861  »  3.8  x  10"6  in 

qm/qNM  *  2.7  x  10“7  (i.e.,  conduct  ng  surface  assumption 

Justified) 

m  -  ~h---A-1 -  1.14  x  10~9  in3/in 
H  ♦  c  lTm-T0J 

or  3*7  x  10T®  in  layer  removed  when  sliding  S  feet. 


SECTION  VIII 


CONCLUSIONS 

Presently  accepted  sliding  friction  theories  are  based  upon  the 
concept  chat  microcontacts  under  high  loads  instantaneously  "weld",  and 
friction  is  the  force  to  shear  these  Junctions,  The  subject  of  sliding 
r’ooity  does  not  enter  the  oodel,  and  only  a  few  studies  show  velocity 
e-'U'Ots  on  friction  data.  The  coefficient  of  friction  between  sliding 
pairs  is  simply  given  as  the  shear  stress  of  the  "welded"  Junctions 
divided  by  the  flow  pressure  of  the  softer  surface.  This  number  agrees 
with  friction  cata  only  qualitatively  (surface  contamination, 
oxidation,  etc.).  However,  hign  sliding  speeds  produces  sufficient 
surface  heat  to  melt  the  one  surface  and  form  (in  10-®  seconds)  a 
melted  fluid  layer  between  the  surfaces.  Fluid  mechanics  enables  one 
to  obtain  more  accurate  expressions  for  sliding  friction  and  also 
expressions  which  are  related  to  sliding  velocity. , 

The  friction  theory,  was  first  developed  for  metal  rotating  bands. 
Essentially  the  theory  separates  the  heat,  generated  from  fluid  shear 
stresses,  into  the  melting  of  one  surface  and  into'  conduction  for  the 
non-melting  surface.  It  was  found  that  a  negligible  fraction  of  the 
heat  went  to  melting  (of  the  rotating  band).  Since  conduction  to  the 
barrel  predominates,  it  is  the  surface  thermal  properties  of  thermal 
conductivity  and  diffuslvlty  that  are  Important.  The  significant 
parameter  for  the  melting  band  is  the  melting  point  temperature. 

During. multiple  firing  conditions  it  is  possible  that  hot  propellant 
gases  cause  the  barrel  surface  temperature  to  exceed  the  rotating  band 
melting, temperature.  Then  there  would  be  conduction  from  the  barrel  to 
the  melt  which  would  add  to  the  shear  stress  generating  heat  to  melt 
the  rotating  band.  The  equations  must  be  altered  accordingly. 

The  metal  melt  was  assumed  to  be  Newtonian  in  behavior.  Tentative 
agreement  with  the  data  of  Montgomery  for  copper  rotating  bands  in  gun 
teats  support  this  assumption.  To  further  Justify  the  overall 
approach,  the  theory  was  extended  to  pin-on-disk  devices.  A  comparison 
was  made  between  the  theory  and  the  Franklin  Institute  data  for  high 
speed  sliding  between  several  pairs  at  various  loads  and  speeds.  Good 
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agreement  was  achieved  in  this  more  strenuous  test  of  the  theory.  At 
lighter  loads  and  slower  sliding  speeds  the  pin  coefficient  of  friction 
changed  considerably,  it  wis  found  that  under  these  conditions  the 
melt  thickness  was  on  the  same  order  as  the  surface  root-mean-square 
roughness  and  full  film  lubrication  was  not  achieved,  typical  me.’ t 
thicknesses  are  approximately  10"5  inches.  Although  this  was  not  a 
problem  for  the  very  smooth,  engraved  rotating  bands,,  future  work  on 
melt  lubrication  should  address  the  melt  of  alcrocontacts. 

Non-Newtonian  melt  rheology  substantially  altered  the  theory  when 
applied  to  plastic  rotating  bands.  Generally, film  thicknesses  did  not 
differ  much  from  the  metal  melts,  but  this  was  caused  by  the  limiting 
shear  stress  condition  of  the  high  loads  and  speeds;  Newtonian  theory 
with  the  high  viscosity  polymer  melts  produced  unreasonable  film 
thicknesses.  Again, most  shear  generated  heat  goes  to  conduction.  The 
most  important  parameter  is  the  limiting  shear  stress  of  the  melt  which 
directly  determines  the  coefficient  of  friction.  Viscosity  and  shear 
modulus  enter  into  film  thickness  calculations.  These  three  parameters 
depend  on  the  pressure  and  temperature  of  the  melt.  Unfortunately, this 
data  for  limiting  shear  stress  and  modulus  Is  practically  non-existent 
for  polymer  melts. 
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